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B.  TECHNICAL  SECTION 
I.  Abstract 


^  The  overall  objective  of  this  proposal  is  to  investi¬ 
gate  the  robustness  to  departures  from  independence  of  methods 
currently  in  use  in  reliability  studies  when  coznpeting  failure 
modes  or  competing  causes  of  failure  associ&ted  with'  a  single 
mode  are  present'  in  a  series  system.  The  first  specific  aim 
is  to  examine  the'^error  one  makes  in  modeling  a  series  system 
by  a  model  whicii  assumes  statistically  independent  conqponent 
lifetimes  whra  in  fact  the  component  lifetimes  follow  some 
multivariate  distribution'.  -  The  second  specific  aim  is  to  assess 
the  effects  of  the  independence  atssuaiption  m  the  error  in 
estimating  component  parameters  from  life  tests  on  series 
S2fstems.  In  both  cases,  estimates  of  such  errors  will  be  deter¬ 
mined  via  mathematicaJ.  analysis  amd  coaoputer  simulations  for 
severail  prominent  multivauriate  distributions.  A  graphical  dis¬ 
play  of  the  errors  for  representative  distributions  will  be 
maide  available  to  researches  who  wish  to  assess  the  possible 
erroneous  assus^tion  of  independent  con^eting  risks.  A  third 
aim  is  to  tighten  the  bounds  on  estimates  of  conponent  reliar- 
bility  when  the  risks  belong  to  a  general  d^endence  class  of 
distributions  (for  example',  positive  quadrant  dependence,  ]^si- 
tive  regression  dependence,  etc.)..  Major  decisions  involving 
reliability  studies,  based  on  competing  risk  methodolpgy,  have 
been  made  in  the  past  and  will  continue  to  be  made  in  the  future. 
This  study  will  provide  that  user  of  |  such  techniques  with  a 
clearer  understemding  of  the  robustness  of  the  analyses  to  de¬ 
partures  from  indepmdent  risks,  an  iassunption  commonly  made 
by  the  methods  currently  in  usS. 
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ZZ.  Specific  Objectives 

The  overall  objective  is  to  investigate  the  robust¬ 
ness  to  depeurtxires  from  independence  to  methods  currently  in 
use  in  reli2d>ility  studies  when  competing  failure  modes  or  com¬ 
peting  causes  of  failure  2U5Sociated  with  a  single  mode  are 
present  in  a  series  system.  We  shall  also  refer  to  such  com¬ 
petitive  events  as  competing  risks.  The  approach  will  be  through 
the  investigation  of  cert2dn  aspects  of  specific  parametric  multi¬ 
variate  distributions  or  by  classes  of  distributions  which  are 
appropriate  in  reli2Ubility  analyses  when  there  are  competing 
risks  present. 

The  specific  objectives  are: 

% 

1}  to  assess  the  error  incurred  in  modeling  system 
life  in  a  series  system  assumed  to  have  indepen¬ 
dent  con^nent  lifetimes  when  in  fact  the  com¬ 
ponent  lifetimes  are  dependent. 

2}  to  assess  the  error  in  estimating  component  param¬ 
eters  (i.e. ,  component  reliability,  mean  com¬ 
ponent  life,  etc.)  in  a  series  system  employing 
either  parametric  or  nonpar ametric  models  which 
assume  independent  component  failure  tines  when 
in  fact  the  lifetimes  ax&  dependent  and  follow 
some  plausible  multivariate  distribution.*  . 

3)  to  derive  bounds  on  component  reliability  when 
the  failure  modes  axe  dependent  and  fall  in  a 
particuleu:  dependence  class  (e.g.,  positive  quad¬ 
rant  dependence,  positive  regression  dependence, 
etc. } • 

4}  to  develop  tests  of  independence,  based  on  data 
collected  from  series  systems,  by  medcing  some 
restrictive  assumption  about  the  structure  of  the 
systems.** 

*  A  plausible  parametric  multivariate  distribution  will  be 
one  that  satisfies  one  of  the  following  conditions: 

i)  the  distribution  of  the  minimum  of  the  component 
failure  times  closely  approximates  widely  accept¬ 
ed  families  of  system  life  distributions. 

or  ii)  the  marginal  distributions  closely  approximate 
the  distributions  of  component  failure  times  in 
the  absence  of  other  failure  modes. 

**This  objective  has  been  added  to  the  original  objectives  be¬ 
cause  it  answers  a  natural  question  raised  by  our  preliminairy 
investigation. 
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ZZX.  '  Introdootion  to  Problen  and  Slgnificaince  of  Stndy 

•  • 

•  •  •  •  • 

Alvin  Weiabe^  (1978)  is  an  editorle^  connent  in  t^' 

pobliahaid  psooeadings*  of  a  voadcahc^’  on  Environmental  l^ologi- 

cal  Hazards  and  Coopating  noted  tbat  ’'tbef  question  of 

****"P*»*"g  risks  will  not  quiefclz^  go  aways  eozsctetions  fox  .eomr^ 

risks  gbcwlfl  be  applied  zmtinely  to  data.*  The  problem 

of  cospeting .  risks  commonly  arises  is  .a  wide  rangd  of  experi- 

mental  situations.  Although  we  shall  confine  our  attention 

in  the  following  discussion  to  those  sitnationa  involving 

series  systans'  in  which  conpeting  failure  modes  or  coopeting 

causes  d  failure'  associated  with  a  single  mode  are  present# 

it  is  certadsly  true  that  we  ‘might  just  as  easily  speak  of 

trials#  animal  experiments#  or  other  medieal' and  bio* 

'logical  studies  where  competing  events  internet  our  study  of 

the'  'main  event  of  interest  (ef  •  Lagokos  (1979))  . 

•  •  • .  •  •  ,  * 

consider  electronic  or  me^anical  systems#  such  as 
satellite  traasaisaion  equipment#,  computers#  aircraft#  missiles 
and  otfasir  ve^pcasy  consisting  of  .several  cooponents  in  series. 
Usnally  each  component  will'  have  a  random  life  length  and  the 
life  'of  the  entire  system  will  end  with  the  failure  of  the 
shortest,  lived  ccmponent^  We 'will  examine  two  situations  more 
cloeely  in  which  eoapeting  risks  play  a  vital  role. 

First#  suppose  we  are 'attempting  to  evaluate  's3rstesi  life 
from  knowle^e  of  the  individual  conponent  lifetimes.  Such  , 
an  avalnatim  will  utilize  either  an  an30.ysis  ixcvolving  math^ 
atatisties  or  a  conputar  simulation. ..  At  a  reccuit 
canfetenba-.'on  Modeling  and  Sisdlation#  HCLean  (1981)  presented 
a  schiiBe  to  sianlate*  the'  life  of  a  •missile  which  consisted  of 
BUizqf  aiajor  eoaponents  in  series^  *  The’  failure  distri.bntion  asso* 
elated  with  eaeb  'eanpeaent  was  usumed  to  be  loaown  Cusually 

or  Weibull.}  *.=  To  aradLve  at  the  systenf  failure  dis* 
tziJbntion#  the  eosponents  wsoce.  assumed  to  act  indep'endehtly  of 
each  'other.  Bealistieally#  this  may  or  may  not  be  the  ease.  , 

Zf  the  'eomponeat  lifetimes  w«e  dspeadent  for  any  reason#  the 
tu-migmitmA  system  failure  distribution  (as  well  as  its  subsequent 
pirimeren  such  'as  systasi  mean  life  and  system  reliability  'for 
a  specified  time)  would  only  crudely  mpproodmate  the  true  . 
dlstributian.  The  first  specifie  aim  of  this  proposal  is  to 
ascertain  the  Incurr^  in  modeling  system  life  in  a 

ffSTi  w  system  assumed  to  have  indapendent  conponent  lifetimes 
(i.e.#  risks)  when#  in  fact#  the  risks  are  dependent. 

Second#  suppose  we  wish  to  evaluate  some  aspect  of  the 
distribution  of  a  particular  failure  mode  based  on  a  Q^ical 
life  test  of  a  series  systeau  The  response  of  interest  is  the 
until  failure  of  a  pairtieular.  mode  of  interest.  Frequently 
response  will  hot  be  observable  due  to  the  occurrence  of 
some  'other  event  which  precludes  failure  associated  with  the 
mode  of  interest.  We  shall  tazm  such  eoapeti^  events  which 
interrupt  our  study  of  the  failure  modes  of  interest  as 

cospeting  risks. 
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Cosopetiag  risks  arise  in  sneh  reliability  studies  when 

1)  the  study  is  temlnated  due  to  a  lack  o£  funds  or  the 
pre^eterBined  period  of  observation  has  eaqpired 
(Type  Z  censoring)  • 

2)  the  study  is  terminated  due  to  a  pre-detexained  number 
of  failnres  of  the  particular  failure  mode  of  interest 
being  observed  (Type  ZI  censoring) . 

3)  some  systssis  fail  because  components  other  than  the 
one  of  interest  malfunction, 

4)  the  coagponent  of*  interest  fails  from  some  cause  other 
than  the  one  of  interest. 

Zn  all  four  situations,  one  may  think  of  the  main  event  of 
interest  as  being  censored,  i.e.,  not  fully  observable.  Zn  the 
first  two  situations,  the  time  to  occurrence  of  the  event  of 
interest  should  be  indepead«it  of  the-  censoring  mechaikisa.  Zn 
su^  Instances,  the  methodology  for  estimating  relevant,  reliabili¬ 
ty  probabilities  has  received  considerable  attention  (cf .  David 
and  Moeschberger  (1978)*,  Kalbfeish  and  Prentice  (1980) ,  Elaztdt- 
johnson  and  Johnson  (1980) ,  Mann,  Schafer,  Singpurwalla  (1974) 
and  Barlow ' and  Proschan  (197S)  for  refercmces  and  discussion)*- 
In  the  third  situation,  l^e  time  to  failure  of  the  coaponent  of 
interest  may  or  may  not  be  independent  of  the  failure  times  of 
’Other  codBqponents  in  the  system.  For  example,  thcure  may  be 
common  environmental  factors  sneh  .as  extreme  temperature  ^diich 
may -affect  the  lifetime  of  several  compoxieats.  Thus  the  question 
of  dependent*  -cospeting  risks  is  raised.  A  •similar  observation 
may  be  made  with  respect  -to  thc^  fourth  situation,  viz.,  failure 
tjawts  associated  with  different  fidJnre  modes  of  a  single  com¬ 
ponent  may  be  dependent.  For  a  very  speci^d.  type  of  dependence,, 
the  models  discussed  by  Marshall-Olkin  (1967),  Lan^berg,  Proschan 
and  Quinsy  (2378) ,  snd  Langberg,  Proschan,  and  Quinsy  (1981) 
allow  one  .to'  convert,  dependent  models  into  independent  ones. 

Zf  no  assusptians  whatever  are  made  about  the  type  of  * 
dependence  between  the  distribution  of  potential  f^uLlure  times, 
there  appears  to  be  little  hope  of  estimating  relevant  eoa^oient 
parasMters.  Zn  soma  situatioas,  one  may  be  i^preciably  misled 
(of.  Tsiatis  (1975),  Peterson  (1976))*  However,  as  Easterling 
(1980)  so  clearly  points  out  in  his  review  of  Birnbaum's  (1979) 
monograph 

"there  seems  to  be  a  need  for  same  robustness 
studies.  How  far  might  one  be  off,  quantita- 
-fcively,  if  his  analysis  is  based  on  incorrect 
assumptions?” 

The  second  specific  aim  will  address  this  important 
issue.  First  if  a  specific  parametric  model  which  assumes 
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independn-t  risks  has  basn  used  in  the  analysis,  it  would 
be  of  interest  to  know  how  the  error  in  estimation;  is 
affected  by  this  assm^tion  of  independence.  That  is,  if 
ijBdependent  specific  parametric  distributions  are  assumed 
for  tixe  failure  times  associated  %rith  different  failure 
modes  when  we  really  should  use  a  bivariate  (or  multivariate 
distribution,  then  eiat  is  the  .magnitude  of  the.  error  in 
estimating  •conponent  parameters?  Secondly,  one  aiay.wish  to 
allow  for  a  less  stringent  type  of  model  assua^ption,  and  ask  ^ 
the  saxDie  question  with  regard  to  the  estimation  error.  That  * 
is,  if  a  nooparametric  analysis  is  performed,  aissuming  in¬ 
dependent  risks,  wfattk  some  types,  of  dependencies  may  be 
psDesent,  what  is  the  magnitude  of  the  estimation  error? 

The  third  specific  ecim  will  attempt  to  (detain  bounds  on 
the  component  reliability  when  the  failure  times  belong  to 
a  broad  dependence  class  <e.g.,  association,  positive  qaadrant 
dependence,  positive  regression  dependence,,  etc.}.  More 
details  will  be  presented  in  the  methods  section.' 

Zn  summary, '  competing  risk  analyses  have  2>een  pmrfoausd 
in  the  past  will  continue  to  be  performed  in  the  ^tnxe. 
Thlg  study  will  provide  the  user  of  such  techniques  with  a 

4  wg  of  the  robustness  to  departures  firom  . 
independent  risks,  .an  assunqition  which  most  of  the  methods 
currently  in  use  assume. 


i 


a":'" 


Moeschberger ,  Melvin  L.  8. 


t 


■•V/ 
;  * » 


i 


6 


$ 


.y, 


i:: 


}t 


•Hf 


IV.  Progress  Report  on  Third  Year's  Work 

A  summary  of  the  first  and  second  year's  work  is  reported  in  the  annual  reports 
dated  October  26th  1983  and  October  26th  1984  respectively.  We  believe  that  during  the 
past  three  years  we  hove  mode  substantial  progress  in  dealing  with  the  objectives  as 
outlined  on  page  4.  In  addition  to  the  papers  and  articles  referred  to  in  the  first  two 
annual  reports,  we  would  like  to  mention  some  of  the  more  recent  work. 

First,  the  recently  published  paper  which  investigates  the  problem  of  improving  the 
product-limit  estimator  of  Kaplan  and  Meier  (1958)  when  there  is  extreme  independent 
right  censoring  is  presented  in  Appendix  A.  This  paper  looks  at  several  techniques  for 
completing  the  product  limit  estimator  by  estimating  the  tail  probability  of  the  survival 
curve  beyond  the  largest  observed  death  time.  Two  methods  are  found  to  work  well  for  a 
variety  of  underlying  distributions.  The  first  method  replaces  those  censored 
observations  larger  than  the  biggest  death  time  by  the  expected  order  statistics, 
conditional  on  the  largest  death,  computed  from  a  Weibull  distribution.  The  Weibull  is 
chosen  since  it  is  known  to  be  a  reasonable  model  for  survival  in  many  situations. 
Parameters  of  the  model  are  estimated  in  several  ways,  but  the  method  of  maximum 
likelihood  seems  to  provide  the  best  results.  The  second  method  replaces  the  constant 
value  of  the  product  limit  estimator  beyond  the  last  death  time  by  the  tail  of  a  Weibull 
survival  function.  Again  parameters  are  estimated  by  a  variety  of  methods  with  the 
maximum  likelihood  estimators  performing  the  best. 

Second,  a  paper  which  obtains  bounds  on  the  component  reliability,  based  on  data 
from  a  series  system,  for  the  Oakes  (1982)  model  has  been  revised.  Since  this  model  has 
the  same  dependence  structure  as  the  random  effects  model  with  w  having  a  gamma 
distribution,  these  bounds  are  good  for  a  general  class  of  distributions.  The  bounds, 
which  are  determined  by  specification  of  a  range  of  coefficients  of  concordance,  are 
found  by  solving  a  differential  equation  in  the  observable  system  reliability  and  crude  life 
on  one  hand  and  the  unobservable  component  survival  function  on  the  other  hand.  This 
revision  is  reproduced  in  Appendix  B. 
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Third,  we  have  submitted  an  overview  paper  for  publication  which  summarizes 
some  of  the  work  performed  during  the  past  three  years.  The  results  of  this  paper  were 
presented  in  an  invited  talk  to  the  Eastern  North  American  Region  of  the  Biometrics 
Society  at  Raleigh,  North  Carolina  in  the  Spring  of  1985.  (See  Appendix  C  for  a  copy  of 
this  paper.) 

Finally,  a  paper  has  been  developed  which  discusses  some  general  properties  of  a 
random  environmental  stress  model.  Estimation  of  parameters  under  the  Gamma  stress 
model  is  considered,  and  a  new  estimator  based  on  the  scaled  total  time  on  test 
transform  is  presented  (See  Appendix  D).  Part  of  the  results  In  this  paper  were  presented 
at  the  International  Statistical  Institute  meeting  in  Amsterdam  in  August,  1985.  A  copy 
of  that  contributed  paper  is  found  in  Appendix  E. 

V.  Methods 

We  refer  to  pages  8-52  of  the  original  proposal  for  a  discussion  of  the  general 
methodology. 
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SuMMAav 

When  there  is  extreme  oenaoring  00  the  r^it,  the  Xapiaa-Meicr  product^imit  estimator  is  known  to 
be  a  biased  estimator  of  the  survival  fluictioii.  Several  modificationt  of  the  Kaplan-Meier  estimator 
are  examined  and  compared  with  respect  to  bias  and  mean  squared  error. 


1.  latradaetiM 

In  human  and  animal  survival  studies,  as  well  as  in  life-testing  experiments  in  the  physical 
sciences,  one  method  of  estimating  the  underlying  survival  distribution  (or  the  idiahility 
of  a  piece  of  equipment)  which  has  received  widespread  attention  is  the  Kaplan-Meier 
prodt^-limit  estimator  (Kaplan  and  Meier,  19S8). 

For  the  situation  in  which  the  longest  time  an  individual  is  in  a  study  (or  on  test)  is  not 
a  feilure  time,  but  rather  a  censored  observation,  it  is  well  known  that  there  are  many 
comi^ex  problems  associated  with  any  statistical  analysis  (Lagakos,  1979).  In  particular, 
the  Kapl^Meier  product-limit  estimator  is  biased  on  the  low  side  (G^  and  Clark, 
197S).  In  the  case  of  many  censored  observations  larger  than  the  largest  observed  feilure 
time,  this  bias  tends  to  be  worse.  Estimated  mean  survival  time  and  selected  percentiles,  as 
wefl  as  other  quantities  dependem  on  knowledge  of  the  tail  of  the  survival  function,  will 
also  exhibit  su^  biases. 

A  practical  situation  which  motivates  this  study  is  a  large-scale  animal  experiment 
conducted  at  the  National  Center  for  Toxicological  Research  (NCTR),  in  which  mice  were 
fed  a  particular  dose  of  a  carcinogen.  The  goal  of  the  experiment  was  to  assess  the  effects 
of  the  carcinogen  on  survival  and  on  age-specific  tumor  incidence.  Toward  this  end,  mice 
were  randomly  divided  into  three  groups  and  followed  until  death  or  until  a  prespecified 
group  censoring  tune  (280, 420,  or  S60  days)  was  reached,  at  which  time  aU  those  still  alive 
in  a  given  group  were  sacrificed.  Often  there  were  numy  surviving  mice  in  all  three  groups 
at  the  sacrifice  times. 

In  general,  we  consider  an  experiment  in  which  n  individuals  are  under  study  and 
censoring  is  permitted.  Let  («),  • ...  Ami  denote  the  m  ordered  feilure  times  of  those  m 
individuals  whose  failure  times  are  actually  observed  (Ad  <  •  •  •  «  The  remaining 
n  -  m  individuals  have  been  censored  at  various  points  in  time.  It  win  be  useful  to  introduce 
the  notation  Sj  to  denote  the  number  of  survivors  just  prior  to  time  Ayi;  that  is,  .Sy  is  the 
number  of  individuals  still  under  observation  at  time  A;),  includmg  the  one  that  died  at 
A;).  Then  the  Kaplan-Meier  product-limit  estimator  (assuming  no  ties  among  the  Ayi)  of 


Key  words:  Adjusted  Kaplan-Meier  survival  estimation;  Bias  of  survival  Auction;  Life-testing;  Right 
censoring;  Survival  analysis. 
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the  underiying  survival  function,  Pit)  *  Pi<r>  /),  is 

1  for  /  <  1(1) 

AO  *  “  n  (*5/  ”  O/'S'  for  t(j)  <  /  <  /(y+i)  (1) 

/>ml 

.0  for  r  >  Wi) 

for7  «  1, . . . ,  m,  where  t(m*i)  *  if  the  longest  time  an  individual  is  on  study  is  a  censoring 
time  or  Wi)  *  « if  the  longest  time  an  individual  is  on  study  is  a  death. 

This  paper  first  proposes,  in  §2,  some  methods  of  “completing’’  the  Kaplan-Meier 
estimator  of  the  survi^^  function  by  (i)  replacing  those  centred  (rijservations  that  are 
larger  than  the  last  observed  fiulure  time  their  expected  order  statistics;  (ii)  using  a 
Weibull  distribution  to  estimate  the  tail  probability  Pit),  for  r  >  r,;  and  (iii)  employing  a 
method  suggested  by  Brown,  Hollander,  and  Kortw  (BHK)  (1974).  The  second  purpose 
is  to  demonstrate  the  magnitude  of  the  bias  and  mean  squared  error  (MSB)  of  the  Kaplan- 
Meier  estimator  and  to  compare  all  methods  of  “completing’*  Pit)  in  the  context  of  the 
aforementioned  mouse  study,  utilizing  simulated  lifetimes  from  exponential,  Weibull. 
lognormal,  and  bathtub^haped  hazard  function  distributions.  These  results  are  presented 
in  §3. 

2.  Coaipletioa  of  Kaplaa-Meier  Prodact-Llmit  FiSrtmaliBr 

2.1  Expected  Order  Statistics 

One  method  of  attempting  to  “complete”  Pit),  t  >  U,  would  be  to  “estimate”  the  failure 
times  for  those  censored  observations  that  are  larger  t^  the  longest  observed  lifetime.  Let 
rtc  be  the  number  of  censored  observations  larger  than  t(m),  A  theorem  regarding  the 
conditional  distributions  of  order  statistics  states  that  for  a  random  sample  of  size  n  from 
a  continuous  parent,  the  conditional  distribution  of  T(k),  given  u>n-~  n^, 

is  just  the  distribution  oftbe  (u  •  n  nrjth  order  statistic  in  a  sample  of  size  rit  drawn  from 
the  parem  distribution  truncated  on  the  left  at  r  •  (see  David,  198 1,  p.  20). 

For  computational  purposes,  take  tr  as  an  estimate  of  the  (n  nr)th  order  statistic.  Then 
find  the  expected  value  of  the  rtr  order  statistics  from  the  parent  distribution  truncated  on 
the  left  at  tr.  Since  the  Weibull  distribution  with  survival  function  Pit)  •  expi-t‘‘/B)  has 
been  widely  accepted  as  providing  a  satisfactory  fit  for  lifetime  data,  it  seems  reasonable  to 
employ  the  results  of  Weibull  distribution  theory  to  complete  Pit),  t  >  u.  (It  should  be 
not^  that  any  distribution  which  is  reasonable  for  the  specific  situation  may  be  used.)  The 
expected  values  of  Weibull  order  statistics  up  to  sam^e  size  40  for  location  parameter 
equal  to  1  and  shape  parameter  equal  to  .5  (0.5)4<1)8  may  be  found  in  Harter  (1969).  For 
la^  sample  sizes,  he  states  a  recurrence  relation  which  may  be  used. 

To  compute  expected  values  of  the  n,  order  statistics  in  question,  values  for  A:  and  9  must 
be  chosen.  One  apiwoech  is  to  use  the  maximum  likelihood  estimators,  ^  and  6,  computed 
by  using  all  obso^ons  to  estimate  k  and  9.  A  second  approach,  due  to  White  (1969), 
uses  least  squares  estimates  of  k  and  9  obtained  by  fitting  the  model 

-  iUk)  In  +  (l/k)  ln(/f(A>,)l  (2) 

to  the  i<;)’s,  where  Hit(j))  is  the  estimated  cumulative  hazard  rate  at  kj)  obtained  from  the 
Kaplan-Meier  estimator.  In  our  Monte  Carlo  study,  we  found  the  maximum  likelihood 
estimators  performed  better  than  the  least  squares  estimators  in  all  cases.  Consequently, 
the  method  of  least  squares  will  be  dropped  fiom  future  discussion  in  this  paper. 

The  survival  function  for  a  Weibull  random  variable,  truncated  on  the  left  at  te,  is 

Piit)  ■  exp{-(r*  -  ti)/9],  t  >  tc. 


(3) 
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So,  by  the  tbeofcm  oo  onkr  stttisdcs  stated  at  the  beginning  of  this  section,  the  conditional 
distribution  of  T^,  given  T(^j  -  (ii  -  «  -  n,  +  1, . . . ,  n)  will  be  approximated 
by  the  (u  n  -f  ne)th  older  statistic  in  a  sample  of  Rc  drawn  from  (3).  For  simplicity,  let 

>  •  u  -  n  fir,  so  that ;  •  1 . ftr.  Now  the  opected  value  of  the/ch  order  statistic  from 

(3)  is 

•  /)  JT  <  ^  dy  (4) 

where  P(y)^  exp(—y^/t'U  y^{t‘‘~  *  0  and  7};j^  is  the/th  order  statistic  in  a  sample 

of  size  fir.  Equation  (4)  can  also  be  written  as 

E(Tj:^)  -  f«,(^  J  J["  (0z*  +  dz  (5) 

where  P(2)  •  exp(-z^  z  •  ( y/tf ) ‘^*  »  0.  Now  ECTj-.^,)  may  be  crudely  estimated  by 

|«(E(Zy:01*  +  tf|''^  (6) 

where  E(Z>;,^)  is  the  expected  value  of  the  /th  order  statistic  from  a  sample  of  size  fir 
determined  fium  Harter’s  (1969)  tables  or  recurrence  rdation,  and  S  and  ^  are  maximum 
liketihood  estimators  of  i  and  k,  respectively. 

These  fir  expected  order  statistics  may  then  be  treated  as  “observed”  lifetimes 

in  adjusting  (or  “completing”)  the  estimated  survival  function  computed  in  (1).  The  area 
under  the  estimated  survival  function  up  to  tc  remains  unchanged.  The  area  under  the 
extended  estimated  survival  function  baa^  on  the  fit  estimated  expected  order  statistics  is 
then  added  to  the  initial  area  to  obtain  a  more  precise  estimate  of  P(t)  [estimated  order 
statistic  (EOS)  extension]. 

2.2  WeUmU  Maximum  Likelihood  Techniques 

A  straightforward  approach  to  completing  P{t)  is  to  set 

P{t)  ■  exp(-r*/fl)  for  t  >  te.  (7) 

Estimates  of  k  and  0  based  on  all  observations  can  be  obtained  by  either  the  maximum 
likdihood  (WTAIL)  or  the  least  squares  method.  However,  our  study  found  the  completion 
using  maximum  likelihood  estimators  was  always  better  in  terms  of  bias  and  mean  squared 
error. 

Ooe  suggestion  for  ostensibly  improving  this  estimator  would  be  to  “tie”  the  estimated 
tail  to  the  product-limit  estimator  at  tc  Two  methods  were  attempted  to  accomplish  this 
goal  Hrst,  the  likdihood  was  maximized  with  respect  to  k  and  0  subject  to  the  constraint 
that  exp(-r^/0)  >■  ^fc).  This  method  will  be  referred  to  as  the  restricted  MLE  tail  probability 
estimate  (RWTAE.  extension).  Second,  a  scale-shift  was  performed  on  the  tail  probability 
in  (7)  to  tie  it  to  the  product-limit  estimator.  This  method  led  to  higher  biases  and  mean 
squared  errors  of  the  survival  function  and  will  be  dropped  from  further  discussion  in  this 
paper. 

2.3  BHK-Type  Methods 

The  Brown-HoUander^Korwar  completion  of  the  product-limit  estimator  sets 
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where  9*  satisfies  P{tt)  >  exp(~r</6*).  In  the  BHK  spirit  we  tried  to  complete  P(t)  by  a 
WeibuU  function  which  used  estimates  of  k  and  9,  k*  and  9*,  that  satisfied  the  following 
twordations: 


him)  -  txBirtC)/9*) 


and 


hkm-n)  -  e3tp(-tC-i)/®*)- 

The  latter  method  also  led  to  consistently  poor  performance  and  the  results  will  not  be 
presented. 
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3.  A  CoapariMM  of  tkcVariooi  Methods 

A  simulation  study  of  data  such  as  that  collected  at  NCTR  was  performed.  Three  groups 
of  48  lifetimes  were  simulated  with  all  testing  stoiqnng  at  280,  420,  and  560  days, 
respectively,  for  the  three  groups.  Distributions  with  mean  survival  times  of  400,  500,  and 
600  days  were  used.  The  generated  lifetimes  greater  than  or  equal  to  the  sacrifice  time  for 
each  particular  group  were  considered  as  censored.  The  remaining  set  of  observed  lifetimes, 
along  with  the  number  censored  at  the  three  sacrifice  times,  constituted  a  single  sample. 
Far  each  of  the  distributioiis  studied,  1000  sudi  samples  were  generated.  Weibull  distribu¬ 
tions  with  shape  parameters  .5,  decreasing  &ilure  rate,  1,  constant  failure  rate,  and  4, 
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increasiiig  £ulure  rate,  were  used.  Lognormal  disttibuiioiis,  dilute  rate  cfaanges  fiom 
increasing  to  decreasing,  with  first  two  moments  comparable  to  the  above  Weibull  distri¬ 
butions  with  k  *■  1  and  A;  >■  4,  were  also  used.  Finally,  a  bathtub  hazard  modd  of  Glaser 
(1980),  &ilure  rate  changes  from  decreasing  to  increasing,  was  used.  This  distribution  is  a 
mixture  of  an  exponential  of  parameter  X  with  probability  1  -  p  and  a  with 

parameter  X  and  index  3  with  prdMbility  p.  Mixing  parameters  of  p  «  .1  and  p»  A  were 
used. 

The  bias  and  MSB  for  the  estimation  of  the  tail  probabilities,  Le.,  the  completed  portion 
of  the  product-limit  estimator,  were  calculated  for  each  hypothesized  distribution  and  for 
each  competing  method  of  completion.  Since  these  results  were  extremely  similar  to  those 
found  in  estimating  mean  survival  time,  A  **  /o  P{t)  dt,  we  show  only  the  bias  and  MSB 
of  each  competing  estimator  of  /t  in  Td>le  1.  This  also  allows  us  to  demonstrate  the 
magnitude  of  the  bias  and  MSB  of  the  produa-iimit  estimator  of  m.  The  bias  and  MSB  for 
estimating  the  90th  percentile  are  also  presented  for  the  various  estimation  methods  in 
Table  2.  As  one  would  expect,  the  Kaplan-Meier  (K-M)  estimator  performs  considerably 
more  poorly  than  the  othn  estimation  schemes.  The  BHK  extension  does  very  well  if  the 
underlying  distribution  is  exponential  or  lognormal  with  filist  two  moments  compatible 
with  the  exponentiaL  BHK  d^  reasonably  well  for  the  bathtub^haped  hazard  model,  but 
it  performs  very  poorly  for  the  Weibull  with  increasing  fiulure  rate  and  for  the  lognormal 
with  first  two  moments  compatible  with  the  WeibulL 

The  remaining  three  extensions  (BOS,  WTAIL,  and  RWTAIL)  appear  to  be  somewhat 
comparable.  Each  of  them  is  best  under  certain  circumstances  although  many  times  the 
biases  and  MSEs  are  so  close  to  one  another  that  they  are  essentially  equivalent  Only  the 
EOS  extension  has  the  desiiaUe  property  of  never  being  worst  It  usually  is  competitive 
with  the  method  that  is  best  Ordering  the  extrasions  firom  the  standpoint  of  simplicity, 
from  simplest  to  most  complex,  we  have  BHK,  WTAIL,  RWTAIL,  and  EOS. 

In  summary,  the  Kaplan-Meier  estimator  should  probably  be  extended  in  the  presence 
of  extreme  censoring.  The  choice  of  extension  depends  on  one’s  knowledge  of  the 
distribution  of  lifetimes  under  consideration  and  the  extent  of  computer  fedlities  available. 
If  the  data  follow  an  exponential-type  distribution  or  if  no  computer  fecilities  are  present, 
the  BHK  method  is  the  extension  of  choice  due  to  its  simplicity.  If  the  data  exhibit  a 
nonconstant  failure  rate  and  computer  fedlities  are  available,  then  the  RWTAIL  or  EOS 
extensions  seem  to  be  advisable. 
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I.  INTRODUCTION 

A  common  problem  In  survival  analysis  Is  Co  escimace  Che  marginal 
survival  funccion  of  che  cime,  Xt  uncil  some  evenc  such  as  remission, 
componenc  failure,  or  deach  due  Co  a  specific  cause  occurs.  Ofcen  obser- 
vacion  of  chis  main  evenc  of  incerest  is  Impossible  due  Co  Che  occurrence 
of  a  competing  risk  at  some  time  Y  <  X,  such  as  censoring,  failure  of  a 
dlfferenc  componenc  in  a  series  system,  or  death  from  some  cause  not 
related  to  the  study.  Standard  statistical  methods,  which  assume  these 
competing  risks  are  independent,  estimate  che  marginal  survival  function 
by  che  Product  Limit  Estimator  of  Kaplan  and  Meier  (1958) .  This  estimator 
has  been  shown  to  be  consistent  for  Che  marginal  survival  funccion  by 
Langberg,  Proschan  and  Qulnzl  (1981)  when  Che  risks  follow  a  constant 
sum  model  defined  by  Williams  and  Lagakos  (1977)  .  When  che  risks  are 
not  in  Che  class  of  constant  sum  models,  Che  Produce  Limit  Escimacor 
is  inconsistenc  and,  in  such  cases,  the  investigator  may  be  appreciably 
misled  by  assuming  independence. 


In  the  competing  risks  framework  we  observe  T  ■  minlmum(X,Y}  and 
I  *  X(X  5,  Y)  where  X(’)  denotes  the  indicator  function.  Tsiatis  (1973) 
and  others  have  shown  that  the  pair  (T,I)  provides  insufficient  information 
to  determine  the  joint  distribution  of  X  and  Y.  That  is,  there  exists  both 
an  independent  and  a  dependent  model  for  (X,Y)  which  produces  the  same 
joint  distribution  for  (T,I).  However,  these  "equivalent”  independent 
and  dependent  joint  distributions  may  have  quite  different  marginal 
distributions.  Also,  due  to  this  ident if lability  problem,  there  may  be 
several  dependent  models  with  different  marginal  structures  which  will 
yield  the  same  observable  information,  (T,I),  In  light  of  the  consequences 
of  the  untestable  Independence  assumption  in  using  the  Product  Limit 
estimator  to  estimate  the  marginal  survival  function  of  X,  it  is  important 
to  have  bounds  on  this  function  based  on  the  observable  random  variables 
(T,I)  and  some  assumptions  on  the  joint  behavior  of  X  and  Y. 

Peterson  (1976)  has  obtained  general  bounds  on  the  marginaL  survival 
function  of  X,  S(x),  based  on  the  estimable  Joint  distribution  of  (T,I). 

Let  Qj^(x)  ■  P(T  >  X,  I  ■  1),  and  Q2(x)  ■  P(T  >  x,  I  ■  0  )  be  the  crude 
survival  functions  of  T.  His  bound,  obtained  from  the  limits  on  the  joint 
distribution  of  (X,  Y)  obtained  by  Frechet  (1951),  is 

Q^(x)  +  Q2(x)  <  S(x)  <.Qj(x)  +  Q^(0). 

Since  these  bounds  allow  br  any  dependence  structure,  they  can  be  very  wide 
and  provide  little  useful  information  to  an  investigator. 

Fisher  and  Kanarek  (1974)  have  obtained  tighter  bounds  on  S(x)  in 
terms  of  a  dependence  measure  a.  Their  model  assumes  that  simultaneous 
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to  the  occurrence  of  Y  an  event  occurs  which  either  stretches  or  contracts 
the  remaining  life  of  X  by  an  amount  associated  with  a.  That  is, 

P(X  >  x|Y  -  y  <  x)  «  P(X  >  y  +  a(x-y)|Y  >  y  +  a(x-y)>.  A  large  a,  for 
example,  implies  that  a  small  survival  after  censoring  is  the  same  as  a-times 
as  much  survival  if  censoring  was  not  present.  They  show  that  if  o  is  assumed 
known,  then  the  marginal  survival  function  can  be  estimated  from  the 
observable  information.  Also  these  estimates,  S,j(x),  are  decreasing  in  o. 

For  their  bounds,  the  investigator  specifies  a  range  of  possible  values 
<  o  <  so  that  (x)  ^  S(x)  ^  S  (x) . 

Recently,  Slud  and  Rubenstein  (1983),  have  proposed  general  bounds. 

They  show  Chat  knowledge  of  Che  function 


„  /  X  1.  P(x  <  X  <  X  -t-d 
p  (x)  ■  lim  .  <  Y  <  *  +<S 


along  with  the  observable  information  (T,I)  is  sufficient  to  uniquely 
determine  the  marginal  distribution  of  X.  These  estimates  Sp(t)  are 

decreasing  functions  of  p  for  fixed  x.  Their  bounds  are  obtained  by 
specifying  a  range  of  possible  values  pj^(x)  ^  p(x)  ^  P2(x)  so  that  if 
p(x)  is  the  true  function  §p2(x}  £  S(x)  £  Sp^Cx)* 

In  this  paper  we  obtain  different  bounds  'on  the  marginal  survival 
function  by  assuming  a  particular  dependence  structure  on  X  and  Y.  These 
bounds  are  functions  of  the  observables  (T,I)  and  a  familiar  dependence 
measure,  the  concordance  probability  between  X  and  Y.  In  Section  2  we 
describe  this  model  in  detail.  In  Section  3  we  derive  the  bounds  and  show 
consistency  when  the  dependence  parameter  is  known.  In  section  4  these 
bounds  are  compared  to  those  obtained  by  Peterson,  Fisher  and  Kanarek, 
and  Slud  and  Rubenstein. 
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II.  THE  MODEL 

The  dependence  structure  vie  shall  employ  to  model  the  joint  survival 
was  first  introduced  by  Clayton  (1978)  to  model  association  in  bivariate 
lifetableii  and,  later,  by  Oakes  (1982)  to  model  bivariate  survival  data. 

Let  S(x)  -  P(X  >  x),  R(y)  -  P(Y  >  y),with  S(0)  -  R(0)  -  1,  be  the 
continuous  univariate  survival  functions  of  the  death  and  censoring  times, 
respectively.  For  9^1  define  F(x,y)  ■  P(X  >  x,  Y  >  y)  by 

F(x,y)  -  +  {^>  -1  (2.1) 

This  joint  distribution  has  marginals  S  and  R.  As  9-^1,  then  (2.1)  reduces 
to  the  joint  distribution  with  independent  marginals.  For  9^, F(x,y) 
min(S(x),  R(y))  Che  bivariate  distribution  with  maximal  positive  association 
for  these  marginals.  The  probability  of  concordance  is  9/(6  1)  so  chat 

Kendall's  (1962)  coefficient  of  concordance  is  t  •  (9  -  l)/(9  +  1)  which 
spans  the  range  0  to  1. 

This  model  has  a  nice  physical  interpretation  in  terms  of  the 
functions  A(;C|Y  ■  y)  and  X(x|y  >  y) ,  the  hazard  functions  of  X  given  Y  ■  y 
and  X  given  Y  >  y,  respectively.  From  (2.1)  one  can  show  that 

A(x|Y  -  y)  -  9A(x|y  >  y) 
or 

P(X  >  xIy  -  y)  -{P(X  >  x|y  >  y))?  (2.2) 

For  9  > 1  the  hazard  rate  of  survival  if  censoring  occurs  at  time  y  is 
9  times  the  hazard  race  of  survival  if  censoring  does  not  occur  at 
time  y.  This  implies  Chat  Che  hazard  race  after  censoring  occurs  is 


accelerated  by  a  factor  of  6  over  the  hazard  rate  if  censoring  had  not 
occurred.  Also  when  6*1,  (2.2),  reduces  to  the  condition  required  by 
Williams  and  Lagakos  (1977)  for  a  model  to  be  constant  sum  and  hence  for 
Che  usual  product  limit  estimator  of  S(c)  to  be  consistent  (See  Basu  and 
Klein  (1982)  for  details). 

Oakes  (1982)  also  shows  chat  (2.1)  can  be  obtained  from  Che  following 

random  effects  model.  Let  S*(x)  ■  exp  {-  t g 1  ^  ^ 

similarly  defined.  Let  U  have  a  gamma  distribution  with  density 

1 _  ..1 

g(w)aw®“^  e“”  and  conditional  on  W  -  w  let  X,Y  be  independent  with 

survival  functions  {S^x) and  {R*(y)}''.  Then,  unconditionally,  X,Y  have 
Che  Joint  survival  function  F(x,y)  given  by  (2.1). 

For  fixed  marginals  S  and  R  the  joint  probability  density  function, 
f(x,y),  can  be  shown  to  be  totally  positive  of  order  2  for  all  8^1. 

This  implies  that  (X,Y)  are  positive  quadrant  dependent.  In  particular, 
one  can  show  chat  for  S,R  fixed  the  family  of  distributions 
F  -  {F(x,y):  6  1}  is  increasing  positive  quadrant  dependent  in  6  as 

defined  by  Ahmed,  seal.  (1979). 
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III.  BOUNDS  ON  MARGINAL  SURVIVAL 

Suppose  that  X  and  Y  have  Che  Joint  dlstrlbuCion  (2.1)  and  leC 
T  «  min  (X,Y),  Chen  Che  survival  function  of  T  is 


F(T)  -  ir  9^ 


and  Che  crude  density  function  associated  with  X, 


P(T  <  t,  X  <  Y),  is  given  by 

qi(0  -  lF(c)f, 

S®(t) 

where  s(c)  ■  -dS(t)/dt. 

Now  consider  the  differencial  equaCion 


(3.1) 


(3.2) 


s(c)/S®(t)  -  q^(typ(t)f 


(3.3) 


and  suppose  6  is  known.  Then  che  solution  of  (3.3)  for  S(c)  is 


5^(0 


(3.4) 


c  q,(u) 

^  -Jw 


if  e  -  1, 


The  functions  F(’)  and  q^(‘)  are  directly  estimable  from  che  data  one 
sees  in  a  competing  risks  experiment.  Let  T^,  ....  T^  denote  the  observed 
test  times  of  n  individuals  put  on  test  and  let  i  >  1,  . . . ,  n  be  1  or  0 
according  to  whether  che  T^  was  an  observation  on  X^  or  Y^  respectively. 


MTJi.TV 


mfm 
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Proof : 

For  0  >■  1,  Che  result  follows  by  a  cheorem  of  Langberg,  Proschan  and 
Quinzi  (1981).  Suppose  chac  0  >  1.  Note  chat  Qj^(t)  -►  Qj^(c)  a.s.  and 

^  A 

F  (u)  F(u)  a.s.  by  Che  strong  law  of  large  numbers.  Since  S.(c)  is  a 

t  ^  ° 

dQ,(u) 

continuous  function  of  ^  ^  in  the  support  of  F(u),  it  suffices  to  show 

0  iF(ii)r 

^  dQj^(u)  ^  <lQj^(u) 

^  ft  '*  ^  fl~  *«s. 

0  [F(u)r  0  [F(u)l” 


Now,  after  an  integration  by  parts. 


^  dQ,(u)  Q  (t)  . 

^  •  zr — 0-  '  Qi(«)<i(Te  > 


0  (F(u)]''  [F(c)l''  0 


Ql^^)  ^  ^  1^1 

/  iQi(u)  -  Q,(u))d(i  )  +  ;  Qi(u)d(i  ) 


F"(u) 


tF(c)] 


r(u) 


F^Cu) 


'  -r - - '  (Q,(U)  -  Q,(u)]d(i  ) 

(F  (u)]®  0  ^  ^  F®(u) 


.dQ,(u) 

+  f  — = - 

''  ft  ' 

0  F  (u)” 

By  the  dominated  convergence  cheorem 

U.  ]  )  dQ.Cu) 

-  0  0  ,7^  *•’•• 


(3.7) 


d 
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Hence,  applying  Che  above  results  Co  (3.7),  Che  result  now  follows:  // 

To  obtain  bounds  on  Che  net  survival  function  based  on  data  from  a 

competing  risks  experiment,  we  proceed  as  follows.  First,  note  that  from 

(3.S)  it  is  true  Chat  S.Ct)  is  a  decreasing  function  of  6  for  fixed  c. 

Also,  as  0  -►  1^  we  have  Sg(t)  f  exp  (-  j  F~^(u)dQj^(u)). 

0 

which  provides  an  upper  bound.  Notice  that  this  upper  bound  corresponds 
to  an  assumption  of  independence.  As  0  •  one  can  show  that  S^Ct)  +  F(t) 

which  corresponds  Co  Peterson's  (1976)  lower  bound. 

In  practice  the  above  bounds,  with  0  »  1,  *,  while  shorter  than 
Peterson's  bounds,  may  still  be  quite  wide. 

Tighter  bounds  may  be  obtained  by  an  investigator  specifying 
a  range  of  possible  values  for  0.  If  Che  sample  size  is  sufficiently 
large  and  0j^  1  9  i  02.  then  Sg  (t)  <  S(c)  <  Sg  (t).  Specifying  0j^,  02 

is  equivalent  to  specifying  a  range  of  values  “^2  the 

coefficient  of  concordance  T  since  9  <■  (1  •f  T)/(l  >'T) .  Hence  the  primary 
value  of  Sg(c)  is  in  putting  bounds  on  S(c)  rather  than  on  estimation  of 


S(t). 


IV.  EXAMPLE  AND  COMPARISONS 

To  llLustrate  che  bounds  obtained  in  Che  previous  section,  consider 
Che  fflortaliCy  data  reported  in  Hoel  (1972).  The  data  was  collected  on  a 
group  of  RFM  strain  male  mice  who  were  subjected  to  a  dose  of  300  rads  of 
radiation  at  age  5-6  weeks.  There  were  three  competing  risks,  thymic 
lymphoma,  reticulutum  cell  sarcoma,  and  other  causes  of  death.  For 
Illustrative  purposes  we  consider  reticulum  cell  sarcoma  as  the  risk  of 
interest . 

Table  1  reports  the  value  of  S-(t)  for  concordance  t  “  (0  -  l)/(6  +  1). 

0 

The  value  of  SgCt)  at  T  -  0  corresponds  to  Nelson's  (1972)  hazard  rate 
estimator  assuming  independence.  Peterson's  upper  and  lower  bound 
(t  *  1)  are  also  reported  as  are  Fisher  and  Kanerek's  bounds  and  the  Slud 
and  Rubens tein  bounds  for  several  values  which  reflect  a  positive 
association  between  risks. 

From  Table  1  we  first'  note  chat  Peterson's  bounds  are  very  wide. 
Substantial  improveoienC  is  obtained  if  one  assumes  a  non-negative 
dependence  structure  between  risks  (See  Table  2).  Further  tightening 
of  these  bounds  is  achieved  by  assuming  that  T  is  in  the  range  0  to  .5 
where  the  width  of  the  boundaries  is  at  most  about  50Z  of  that  of  Peterson's 
bounds. 

Substantial  improvement  in  the  general  bounds  is  also  obtained  by 
the  bounds  of  Fisher  and  Kanerek  or  Slud  and  Rubenstein.  The  bounds  of 
Fisher  and  Kanerek  assume  a  specific  censoring  pattern  and  require  a 
specification  of  a  stretching  constant  a.  Without  some  additional  informa¬ 
tion,  such  specification  may  be  impossible.  Slud  and  Rubenstein 's  bounds 
are  for  the  general  dependence  structure.  Their  bounds  require  che 
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specification  of  the  p(c)  function.  This  function  is  a  quantity  which 
is  not  easily  conceptualized  by  investigators  from  either  a  statistical 
or  biological  perspective.  This  makes  it  questionable  whether  reasoable 
upper  and  lower  bounds  on  p(t)  can  be  extracted  from  one's  prior  beliefs. 
The  major  advantage  of  the  bounds  printed  in  this  paper  is  that  they 
require  only  the  specification  of  an  upper  and  lower  concordance,  a 
measure  quite  familiar  to  most  investigators  and  easily  explainable  to 
nonstatisticians. 
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Table  2 

ElELATIVE  SIZE  OF  THE  BOUNDS  ON  NET  SURVIVAL 
FOR  AN  ASSUMED  DEPENDENCE  STRUCTURE 
AS  COMPARED  TO  PETERSON'S  BOUNDS 


0  £  T  £  1 

0  £  T  £  .5 

0  <  T 

.9707 

.0879 

.2674 

.9352 

.2449 

.5931 

.7338 

.5171 

.6787 

.6722 

.5120 

.6298 

.5009 

.4420 

.4870 

.3831 

.3576 

.3797 

.2883 

.2767 

.2833 

.0600 

.0600 

.0600 
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ABSTRACT 

This  article  investigates  the  consequences  of  departures  ftta 
independence  t4ien  the  oanponent  lifetimes  in  a  series  system  are 
exponentially  distrihuted.  Such  departure  are  studied  vhen  the 
joint  distribution  is  assumed  to  follow  either  one  of  ^  three 
Gumbel  bivariate  exponential  models,  the  Dounton  bivariate 
exponential,  or  the  Oakes  bivariate  exponential  model.  TVo 
di^inct  situations  are  considered.  First,  in  theoretical  modeling 
of  series  systems,  when  the  distribution  of  the  cOTponent  lifetimes 
is  assumed,  one  wishes  to  conpute  system  reliability  arxl  mean 
system  life.  Second,  errors  in  parametric  and  nonparamepric 
estimation  of  oomponent  reliability  and  component  mean  life  are 
studied  based  on  life-test  data  collected  on  series  systesns 
the  assumption  of  ixKlependence  is  made  erroneously.  Systems  with 
two  components  are  studied. 

KEY  WORDS:  Competing  risks;  Component  life;  Modeling  series 
systems;  Robustness  studies;  System  reliebility; 

Gumbel  bivariete  exponential;  Dounton  bivariate 
exponential;  Oakes  bivariate  exponential. 


I.  Introduction 

Consider  a  system  oonsisting  of  several  ooB^poneixts  linked  in 
series.  For  sudi  a  system  the  failure  of  any  one  of  the  components 
causes  the  system  to  fail.  Ih  a  biological  or  medical  oontext 
we  can  consider  the  components  to  be  different  lethal  diseases 
and/or  different  reasons  for  removal  from  the  study.  In  a  clini¬ 
cal  trials  Aremeuork  the  prinary  response  of  interest*  death  or 
remission,  and  censoring  can  be  considered  as  components  of  the 
system.  This  general  fomulation  has  been  detailed  in  the  theory 
of  oompeting  risks  (cf.  David  and  Moeschberger  (1978)). 

A  oonmon  assumption  in  such  a  fomulation  is  that  the 
component  lifetimes  are  statistically  independent.  Several 
authors  have  shown  that  based  on  data  from  series  systaos  only, 
this  assunption,  by  itself,  is  not  testable  because  there  is  no 
way  to  distinguish  between  independent  or  depenient  ccmponent 
lifetimes  (see  Basu  (1981),  Basu  and  Klein  (1982),  tHUer  (1977), 
Peterson  (1976),  etc.).  However,  several  authors  (see  Lagakos 
(1979)  p.  1S2  and  Easterling  (1980)  p.  131)  have  pointed  out  the 
need  to  determine,  quantitatively,  how  far  off  one  ai^t  be  if  an 
analysis  is  based  on  an  incorrect  askmption  of  independence. 

To  study  effects  of  erroneously  assuning  independence  we 
shall  assune  that  each  of  the  components  have  exponentially  dis¬ 
tributed  lifetimes  when  tested  separately  and  that  the  property  of 
marginal  exponentiality  will  be  preserved  even  though  seme 
dependence  nmy  be  induced  ^en  the  cemponents  are  linked  in  series. 
The  assiApticn  of  exponentially  distributed  component  lifetimes 
has  been  made  by  Mann  and  Grubbs  (197H) ,  when  finding  confidence 
bounds  on  system  reliability,  Boarctani  and  Kendall  (1970) ,  vhen 
estimating  component  lifetimes  from  system  data,  and  Kiyaaura 
(1982) ,  then  combining  con^onent  and  sytem  data.  (See  Barlcw  and 
Proschan  Q975}  or  Mam,  Schaffer,  and  Singpurwalla  (1974)  for  a 


more  oonplete  review. }  Vie  shall  model  the  dependence  structure  by 
the  tivSae  models  of  Gimibel  0.960),  a  model  proposed  by  Do«citon 
097Q),  and  a  model  described  by  Oakes  0982).  These  models  are 
briefly  described  in  Section  2. 

The  effects  of  a  departure  fircm  tiie  assuaption  of  ind^iendent 
ocmponeiit  lifetimes  will  be  addressed  for  tvo  distinct  situations. 
The  first  situation  arises  in  modeling  the  pa?fbnnance  of  a 
theoretical  series  system  constructed  from  tvo  ccnponents.  Here, 
based  on  testing  each  oonponent  separately  or  on  engineering 
design  principles,  it  is  reasonable  to  assune  that  the  components 
are  exponentially  distriJbuted  with  known  parameter  values.  Baised 
on  this  information,  we  wi^  to  predict  parameters  such  as  the 
mean  life  or  reliability  of  a  series  system  oonstructed  from 
these  oomponents.  In  Section  3  we  describe  how  these  quantities 
are  affected  by  d^artures  from  independence. 

The  second  situation  involves  making  inferences  about 
oonponent  lifetiffle  distributions  from  data  collected  on  series 
systems.  Coomonly,  data  collected  on  such  systems  are  analyzed 
by  assimiing  a  constant-sum  model,  of  which  independence  is  a 
special  case  (compare  Williams  and  Lagakos  (1977)  and  Lagatos  and 
Williams  (1978)).  In  Section  4.1  we  study  the  properties  of  the 
maxinum  likelihood  estimators  of  the  component  mean  life  calcu¬ 
lated  under  an  erroneous  assumptiOTi  of  independent  exponential 
component  lifetimes  as  mentioned  above.  Because  of  the  wide 
spr^d  use  of  the  nonparanetric  estimator  proposed  by  Kaplan- 
Meier  (19S8)  for  the  component  reliability  v«  study  in  Section  4.2 
its  pcoperties,  ihen  the  nmrginal  reliabilities  are  exponential 
and,  independenoe  is  incorrectly.  assuaedV 

2.  The  Models 

Consider  a  two  ccmpcnent  series  system  with  component  life 
lengths  Xj^,  Xj.  Suppose  that  each  Xj^  has  an  exponential  aavival 
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function 


C2.1)  fjCtl  s  POCj  >-  t)  3  expC->^  t).  X.  >  0.  t  >  Q. 

This  assuaption  is  node  on  tiie  basis  of  eetansive  testing  of 
each  oonponent  sqiaxately  or  on  loMwledge  of  tiie  underlying 
mechanisn  of  failure. 

Tb  examine  the  effects  of  a  departure  from  independence  we 
consider  five  bivariate  exponential  models »  each  with  marginals 
equivalent  to  (2.1)  •  The  first  three  models  are  due  to  Gumbel 
(1960);  the  last  two  models  are  due  to  Downton  (1970)  and  Oalces 
(1982). 

2.1  Gtanbel's  ?fodel  A 

For  this  model  the  joint  survival  function  is 


(2.2)  P()€j^  >  Xj^,  Xj  >  Xj)  *  ejq>C-  Xj^Xj^  -  X^Xj  - 

Xj^,  *2  >  0*  Xj^,  Xj  >  0,  0  <  Xj^2  1 

The  correlation  betweeh  )C.  >  X.  is 
XX  ^ 

0  3  -  expa^X^/X-^^)  Ej^(-  X^X^/X^^)  -  1, 

12  ^ 

Utere  E.(z)  =  /*_  exp(-u)  du  is  the  integrated  logarithm 
^  u 

For  this  model  0  varies  from  -  .40365  to  0  as  ^12  decreases  from 
Xj^X^  to  0.  It  is  never  positive.  The  regression  X^  on  is  non- 
leanear  with  ^ 

ECX.  1X2  3  X2)  3  (Xj^  ♦  Xj^2^  -  *  Xj^2^)  - 


2.2  Gunbel’s  Model  B 

For  this  model  the  joint  survival  function  is 
C2.3)  .P(Xj^>  Xj^,  Xj  >  X2)  *  expC-  X^^Xj^  -  X2X2){1  ♦  4p(l-exp(-Xj^Xj^)) 
a-  expC-  X2X2)}.  X^,  Xj  >r0,  ^  ^  0,  -  1/4  <  p  <  1/4. 

The  oorvelationy  Pf  may  be  positive  or  negative.  ®e  regression 


of  X2  on  is  agadji  nonlinear  ultfi 

E0Cj^|X2  a  Xjl  a  (1  ♦  2P  -  •«>  expC-  ^2^1^ 

The  effects  of  a  depaz'Cure  fxxan  indepoxtence  on  modeling  systea 
reliabili^  euid  estimating  oompment  reliabilities  has  been  studied 
in  detail  in  Moeschberger  and  Klein  (1984). 

2.3  Gunbel*s  ^txlel  C 

For  this  model  the  joint  survival  function  is 
P(Xj^  >  Xj^f  X2  >  Xj)  =  exp{“  [(Xj^Xj^)™  +  , 

X^*  Xj  >  0,  m  ^  1,  Xj^,  X2  ^  0« 

The  correlation  is 
0  =  (4  ♦  an)  r (cos  0  sin  9)® 

J  Z - 7771*SSr'  de-l 

0  |_6os  0“  ♦  sin  0“) 

«4iich  varies  frca.  0  to  1.  For  this  model  m  a  1  corresponds  to 
independence  and  as  m 

(2.S)  POCj^  >  Xj^,  Xj  >  Xj)  *  ndiumin  Cexp(-  exp(- 

the  Frecbet  (19S8)  upper  bound  for  these  nargindds. 

2.4  Downton*s  Model 

DoMiton  (1970)  suggests  modeling  bivariate  exponential  systems 
by  a  successive  damage  model.  This  model  assunes  that  in  a  tvo 
component  system  the  times  bet«#een  successive  shocks  on  each 
component  have' ind^endent  exponential  distributions' and  that  the 
nunber  of  shodcs  required  to  cause  each  ooupdiott  to  fail  follows 
a  bivariate  geometric  distribution.  The  joint  probability 
d«isity  function  of  the  compcmdit  life  times  is 


C2.5)  f(Xj^^X2}  “  CX^x^ -t  12X2^  ^1qi 


L-p 


I  ^  Io^2^‘pX^^X2 


*1* 

i 


1 

:V 


vtere  IqC*}  is  the  modified  Bessel  function  of  the  first  kind  of 

order  zero,  and  >  Q,  Xj^,  Xj  ^  0,  0  ^  p  ^  The  ootrela- 

tion  between  X^,  X2  is  p  tiiich  spans  the  interval  CQ,l3»  As  P  1 

the  joint  survival  function  of  approaches  the  upper 

F^rechet  distribution  C2.S).  For  this  model 

ECXj^lXj  =  Xj)  s  Cl-  P/Xj,  ♦  P  XjXj/X^. 

% 

2.5  Oakes'  Model 

Oakes  (1982)  has  proposed  a  model  for  bivariate  survival 
data.  This  model  vas  first  proposed  by  Qayton  (1978)  to  model 
association  in  bivariate  lifetables.  Special  cases  of  Oakes' 
general  ■nryU'i  have  been  suggested  by  LindleyandSingpurvalla  (1985) 
and  Hutchinson  (1981 ) . 

For  this  model  the  joint  survival  probability  is 
(2.7)  P(Xj^>  x^,  X2>  x^)  s  Cexp(X^(9-l)Xj^)  ♦  cxp(X2(8-l)x2)-l]"^^“^^ 
where  Xj^,  X^  ^  ®  i 

For  8  s  1,  Xj^,  X2  are  independent  and  PCXj^  >  Xj^,  Xj  >  ^2)  (2.5) 

as  9  For  this  ntxlel  Kendall's  (1962)  coeffcient  of  concordance 

is  T  =  (9-1)/ (9+1)  vhich  spans  the  range  0  to  1.  The  correlation, 
p,  also  spans  the  rai^ge  0  to  1  and  is  found  numerically. 

This  nxxlel  has  the  following  physical  interpretation.  Let 
r(x^|X2  =  ^2)  and  r(x^|X2  >  X2)  be  the  conditional  failure  rates 
of  X^  given  X2  =  X2  and  X2  >  X2.  Then  r(Xj^(X2  =  X2)  =  9  r(xj_|X2  >^2^ 
The  ncxlel  can  also  be  derived  from  a  random  effects  model. 

This  fonoulation  assuoes  that  i4ien  the  conponents  are  tested 
separately  under  conditions  the  component  survival  functions 

ar«  S^(t)  s  expC-  exp(X^t(9-l))  ♦  IJ,  i  *  1,  2.  and  that  vhen  the 
two  oonponents  are  put  in  a  series  system  in  the  operating 
envirorment  there  is  a  random  factor  W  thich  simultaneously  changes 
each  oonponent  life  distribution  to  S^(t) ..  ^If  W  has  a  gpmoa 

distribution  with  density  function  fCx)  ow  e  then. 


unconditional!/,  the  joint  survived,  function  (2.7)  holds. 
2.6  Fr&het  Bounds 


Frechet  0.9S8)  obteiined  bounds  of  the  joint  survived  functions 
vhich  can  be  obtained  for  set  of  marginal  distributions.  For 
exponential  margixvals  these  are 


— X.X.  ”^2^ 

MAXIMUM  (e  ^  ^  *  e  -  1,  01  <  P(Xj^  >  X2  >  x^)  < 

MINIMJM  (e  e  ^  . 

For  this  set  of  marginals  the  lower  FTechiet  distrilxition  has 
correlation  -.634  and  the  t^Jper  FTechet  distribution  has 
correlation  1.0.  These  are  the  mininal  and  maximal  correlations 


for  exponential  marginals. 


3.  Errors  in  Modeling  System  Life 

Sujqpose  that  based  on  extensive  testing  or  based  on  theore¬ 
tical  considerations  eech  of  the  two  components  in  a  series 
system  is  Jcnowi  to  have  an  exponential  distribution,  (2.1)  with 
narginal  means  1/Aj^,  l/Aj*  respectively.  It  is  of  interest  to 
predict  the  system  reliability  F(t)  *  P(Xj^  >  t,  >  t)  and  the 

system  mean  life  W  "  F(t)  dt.  If  the  investigator  assumes  that 
the  two  components  are  independent  then  the  system  reliability  is 

(3.1)  Fj(t)  =  exp(  -(Aj^  ♦  A^lt)  and  system  mean  life  is  Wj  =  l/(Aj^  +  A 

If  the  ccn?»nents  are  not  independent,  but  in  fact  follow  one 
of  the  (ttxiels  in  Section  2,  then  a  measure  of  the  effects  of 

ind^MndgnOa-Aift)  *  (FCt)  -  Fj(t))/FjCt)  and 
5  r  (|i  -  Hji/Uj, '  for  prodbctinff  sy^aa  reliability' and  system  mean 
life  ^respectively,  »iiere  FCt)  and  ii'are  oonputed  under  the 
appropriate  dependent  model.  Values  of  FCt)  can  be  cooputed 
directly  from.  (2.2),  (2.3),  (2.4),  (2.6)  (by  nuherical  integration) 
or  ftom  (2.7).  Expressions  for  w  are  ^ven  in  Appendix  11  AH' 


expressions  for  A(t}  and  ^depend  on  the  values  of  and  Xj  only 
tiuxxjgh  the  ratio  X2/X2  ^  ^  ^or  K  <  1  tiie  values  are  equiva¬ 
lent  to  those  for  IC*  s  ^  For  the  upper  Frediet  distrihution« 

K 

-  TO" 

ACtp)  a  p  -  1  vhere  K  ^  1  and  tp  ~  the  point 

‘>*'ere  -  P*  Also  £  =  1/K  for  K  >  1.  For  tiie  lower  Fbechet 

distribution 

-  1  K  K  1 

AvtpJ  -  p  +p  — p  -llXp  ♦p-l>0 

-  1  othendse 

and  6  =  •»•  K  ♦  1  -  (K+1)^  Y  ♦  (m)  UCi)  vhere  Y  is  the  solu- 

K  ~ir“ 

tion  of  the  equality  >  X  a  1.  Table  1  gives  the  values  of 
ACtp)  X  100%  and  &  X  100%  for  p  a  .9,  .7,  .5,  .3,  .1  for  the  upper 
and  lotasr  F)?&^et  distributions. 

Frcxa  Table  1  we  see  that  the  largest  percent  error  occurs 
when  the  parameters  are  equal.  Also  for  fixed  K  there  is  rela¬ 
tively  snail  error  in  estimating  system  reliability  by  modeling 
a  dependent  system  by  an  independent  system  vhen  FCt)  is  large* 

For  snaller  values  of  system  reliability  one  can  be  appreciably 
misled.  Errors  in  estimating  system  mean  life  appear  to  be 
substantial  unless  one  component  has  considerably  longer  marginal 
life  than  the  second  one.  In  that  instance,  one  can  see  instinc¬ 
tively  that  the  correlation  would  have  a  mininal  impact. 

Figures  lA-lC  and  2Ar2C  ace  plo^  of  A^tpl^for.p  s  .24,.  .5,  .75 
X^  s  1  and  X2  =  1*  1«S  for  the  five  models. described  in  Sectiim  2. 
Figures  3A,  3B  are  plots  of  £  fbr^all  five  models  as  a..function  of 
the  ooirelation.  Item  these  plot»  note  tiiat  for  positive  correla¬ 
tions  the  most  nodeling  error  occurs  for  the  Guabel  C  model.  For 
relatively  small  caErelation,--.2S.  <,  p  £  .25  there  may  still  .be  a 


modex>ate  nodeli^  error,  on  the  order  at  least  ^  10%  for  pre¬ 
dicting  system  reliability  at  F^Ct)  s  .25,  or  F^Ct)  s  .$  and  fbr 
estinating  the  mean  system  life. 


TABLE  1 

UPfCR  ANB  LQHEft  BOUNDS  ON  THE  PERCENT  ERROR  IN  MSaiNS  STSTEH  LIFE 
F(T)-0.t  F(T)«0.7  FITI^.S  f(T)-«.J  F(TI«0.1  HEAR  LIFE 


loneT" 

BOUND 

UPKR 

aOUNO 

LONER  “ 
BOUND 

UPPER 

BOUND 

LONER 

BOUND 

UPPER 

BOUND 

LONER 

BOUND 

U^  LONER 
BOUND  BOUND 

UPPER 

BOUND 

LONER 

BOUND 

UPPER 

BOUND 

•0,21 

3.41 

-1.81 

11.32 

-17.16 

41.42 

-100.00 

82.37  -100.00  216.23 

•38.63  100.00 

•0.2A 

3.37 

-1.31 

12.62 

-13.27 

23.11 

•100.00 

41.38  -100.00  113.44 

•n.o7 

30.00 

-0.22 

2.A7 

•2.86 

1.33 

-12.10 

18.12 

-31.32 

33.12  -100.00 

77.83 

•34.83 

33.33 

-0.11 

2.13 

-2.44 

7.31 

-11.02 

14.87 

-44.11 

27.23  -100.00 

58.41 

•32.83 

23.00 

•O.U 

1.77 

-2.12 

6.12 

-1.57 

12.23 

'38.38 

22.22  -100.00 

46.78 

•31.06 

20.00 

-0.14 

1.32 

-1.87 

3.23 

-8.43 

10.41 

-33.11 

18.n  -100.00 

38.13 

-21.33 

16.67 

,  -0,13 

1.33 

-1.67 

4.36 

-7.33 

1.05 

•30.33 

16.24  -100.00 

33.33 

•28.18 

14.21 

-0,12 

1.10 

-1,31 

4.04 

•6.82 

8.01 

-27.42 

14.31  -100.00 

21.13 

-26.91 

12.30 

-O.U 

1.06 

-1.17 

3.63 

-4.22 

7.18 

-23.02 

12.71  -100.00 

23.81 

-23.13 

11.11 

-0,10 

0.16 

-1,26 

3.30 

-3.71 

6.30 

-22.11 

11.37  -100.00 

n.28 

-24.98 

10.00 

-0.01 

0.88 

-1.17 

3.02 

-3,28 

3.15 

-21.27 

10.33  -100.00 

21.13 

•24.12 

1.01 

•0.00 

0.81 

-1.08 

2.78 

-4,11 

3.48 

-11.78 

1.70  -100.00 

11.38 

-23.34 

8.33 

-0.00 

0.76 

-1.01 

2.38 

-4,51 

3.08 

.  -18.41 

8.18  -100.00 

17.88 

•22.62 

7.61 

-0.07 

0.70 

-0.13 

2.41 

-4.30 

4.73 

-17.33 

8.36  -100.00 

16.39 

-21.96 

7.14 

-0.07 

0.66 

-0,10 

2.23 

-4.05 

4.43 

-16.35 

7.82  -100.00 

13.48 

-21.33 

6.67 

-0.06 

0.62 

-0.83 

2,12 

-3.83 

4.16 

-13.45 

7.34  -100.00 

14.30 

-20.78 

6.23 

-0.06 

0.31 

-0.80 

2.00 

-3.63 

3.13 

-14.63 

6.92  -100.00 

13.63 

-20.26 

3.88 

-0.06 

0.36 

-0.76 

1.81 

-3.43 

3.72 

-13.13 

6.34  -100.00 

12.88 

-11.76 

3.36 

-0.06 

0.33 

-o.n 

1.80 

-3.21 

3.33 

-13.27 

6.20  -100.00 

12.20 

-11.30 

3.26 

-0.03 

0.30 

-0.61 

1.71 

-3.14 

3.36 

'12.67 

3.10  -12.26 

11.31 

-18.87 

3.00 

-0.03 

0.48 

•0.66 

1.63 

-3.00 

3.20 

-12.13 

3.63  -88.34 

11.03 

-18.4 

4.76 

-0.03 

0.46 

.  -0.64 

1.36 

-2.88 

3.06 

'11.63 

3.37  -84.73 

.10.33 

-18.08 

.  4.33 

-0.03 

0.44 

•0.61 

1.30 

-2.76 

2.13 

•11.16 

3.14  -81.41 

10.07 

-17.71 

4.33 

•0.04 

0.4» 

-0.31 

1.44 

2.81 

-10.74 

4.13  -78.34 

1.63 

-17.37 

4.17 

-0.04 

0.41 

-0.57 

1.38 

•2.36 

2.70 

-10.34 

4.74  -73.41 

1.26 

•17.04 

4.00 

4.  Errors  in  Estimating  Gon^nnent  Paxometers 
4.1  Paiometric  Estirtation 

In  this  section  we  examine  tive  effects  of  incorrectly  assuming 
independence  on  the  magnitude  of  the  estimation  error  in  esti- 
neting  tdie  first  ocmportent  mean  life  based  on  data  from  series 
systems.  Suppose  that  n  series  systems  are  put  cwi  test.  For 
each  systan  we  observe  the  system  failure  time  and  which  component 
caused  the  failure.  Let  n*  dentoe  the  number  of  systems  where 
the  system  failure  was  caused  by  failure  of  the  i^  component, 
i  =  1,  2,  and  let  T  be  the  total  time  on  test  for  edl  n  systems- 
If  we  assune  that  tite  oompor^t  lifetimes  are  independent  and 
exponentially  distributed  then  Moeschberger  and  David  (1971)  show 
that  ■dte  maxinua  likelihood  estisator  of  Uj_»  tiie  first  oomportent 
mean  life  is 

A 

C4.1)  s  for  n^^  >  0. 

This  estimator  is  asymptotically  unbiased  and  for  n  fixvite 

A 

E(uj^)  s  E(T) 'ECl/nj^lnj^  >  0)  due  to  the  independence  of  T  and  n^. 

Suppose  now  that  the  two  component  lifetimes  are  not  inde¬ 
pendent  but  follow  one  of  the  models  discribed  in  Section  2.  If 
we  incorrectly  assume  independence  then  a  measure  of  the  excess 
bias  due  to  incorrectly  assuming  independence  is 
B  =  CECuj^ I  Dependent  model}  -  ECUj^l  independence) 
the  dependent  models  under  consideration  T  and  are  independent. 
For  large  n,  B  converges  to  (y/p  -  where  w  is  the  mean 

system  life  and  p  is  the  probability  the  first  cooqponent  fails 
first,  computed  under  the  ^pendent  model.  For  finite  n, 

E(uj^)  s  n  w  Ep(l/nj^|nj^  >  Q)  ocmputod  under  the  appropriate  model, 

where  EpO/oj^lnj^  >  Q)  =  ^I^(JJ)p^a-p)*^/K  /  U-Cl-p)").  Expressions 
for  u  and  p  are  given  in  Appendix  1  and  Appendix  2,  respectively. 


The  expressions  d^iend  on  A2  only  throu^  tive  ratio  K  =  A^/^2 
Fbr  all  nodels,  p  a  1/2  tiien  K  a  1. 

For  tile  upper  Fr^chet  distribution  p  a  q  if  K  <  1;  1/2  if  K 
and  1  if  K  >  1.  Kance  for  K  <  1  no  failures  from  the  first  ooB^* 
parent  are  ever  observed  so  that  the  modeling  error  B  becomes 
infinite  for  all  n.  For  K  >  1,  p  =  1  and  u  =  so  that 
B  s  Cl  -  I  Independence) /Uj^)  viiich  tends  to  0  as  n  In 

this  case  the  models  with  oorrelation  ranging  from  0  to  1  have  B 
increasijTg  for  p  <  Pg  and  decreasing  for  p  >  pg.  For  the  lower 
Frechet  distribution,  p  is  the  value  of  X  vhich  solves  the 
equation  X  -  1  a  0.  For  K  <  1  wehave  p  <  1/2  and  for  K  >  1 
we  have  p  >  1/2.  Table  2  gives  'Oie  value  of  B  for  n  ~  2S,  SO,  <■ 
for  the  two  FVechet  distributions.  It  ailso  gives  the  maDcinun 
modeling  error  for  the  Gumbel  C  model  «hixdi  is  an  indication  of 
nBxiinal  excess  modeling  error. 


nxm  Table  2  vte  note  that  the  dependence  structure  exerts  a  large 
effect  on  esthnating  the  analler  of  the  two  oomponent  neans  ennd 
that  either  effect  is  most  eseaggerated  fbr  snail  saa(d.e  sizes. 

For  K  ^  1  tiiere  is  vecy  little  sample  size  effect  on  tiie  modeling 
error.  For  K  strictly  bigger  than  one  the  maximun  bias  under  the 
Gunbel  C  model  decreases  with  K  and  the  correlation  at  vliich  this 
maxiflun  is  attained  also  decreases  to  0.  Figures  4A-4C  are  plots 
of  B  as  n  •  for  K  =  3/2,  1,  2/3,  respecitvely.  Figures  SA-SC 
are  plots  of  B  for  n  =  10  for  K  =  3/2,  1,  2/3,  respectively. 
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Figures  SA-SC  for  p  =  .25,  6A-6C  for  p  =  .5  and  7A-7C  for 
p  =  .75  are  plots  of  Aj^(p)  for  the  5  models  and  k  =  3/2,  1,  2/3. 
As  in  the  previous  figures  one  can  see  that  for  even  a  small 
departure  from  independence  the  relative  effect  of  dependence 
can  be  quite  large. 


4.2  Nonparantetric  Estijution 

A  second  approach  to  the  probLen  of  estimating  ocnponent 
paz>ameters  is  via  the  nonpazemetric  estimator  of  lQ^>lan  and 
Meier  C19S8}.  investigators  tte  routinely  use  nonparametric 
techniques  may  take  this  approa^  in  hopes  of  obtaining  estimators 
that  are  robust  with  respect  to  the  assunption  of  exponentiality. 
However  this  estimator  is  not  necessarily  robust  to-  the  assunption 
of  independence. 

The  product  limit  estimator,  assuning  independent  risks  is 
constructed  as  follows.  Suppose  tijat  n  systems  are  put  on  test 
and  let  r^,  ...,  r^^  be  tite  ranks  of  tile  ordered  n^  failures 
from  cause  i,  ^iQ)*^-***  order  lifetimes. 

The  estimator  of  the  oaaponent  reliability  for  the  i^  component  is 

4.2.1  *SiCx).lif  x<Xijy, 

j(i.x)  x>  X(£Q) 

"  "  "  ‘ 

j*l  n  -  rj^.  ♦  1 

jCi,x)  is  the  largest  value  of  j  for  which  <  x.  This 

estimator  is  asynyotically  unbiased  vAmn  tiie  component  lifetimes 
are  independent. 

When  the  risks  are  dependent  Klein  and  Hoeschberger  (1983) 

A 

show  that  S^Ct)  is  not  estimating  tiie  marginal  conponent  relia¬ 
bility,  but  rather  it  is  estimating  consistently  another  survival 
function 

4.2.2  t^Cx)  =  exp^-l"^  d  Q^Ct)  ^  where. 

^Ct)  s  PCminimum  CX^,X2}  >  t)  and  Q^Ct)  s 

PCminOCj^.Xj}  <  t,  minO^iXj)  ^  Xj.),  i  ^  l,  2.  Expressions  for 


F^Ct)  for  tile  five  nodels  of  iatecesi;  ere  given  in  Appeniix  3. 

A  measure  of  the  affect  of  dependence  in  using  the  Product 
limit  estiaator  with,  dependoit  risks  is  d^Cp)  - 
where  t^  is  the  time  ^ere  the  true  ocnponent  reliability  is  p. 
Aj^Cp)  is  again  only  a  function  of  k  »  ^  “PP*r 

Frechet  distribution 

Aj^Cp)  s  p~^  -  1  for  k  <  1 

p“^^^  -  1  for  k  s  1  since 

Q  for  k  >  1 

for  k  <  1  i^Ct)  =  1  for  all  t  since  the  first  component  never 
fails,  while  for  k  >  1  all  failures  are  due  to  the  first  oomponent 
For' those  nodels  with  ootrelation  spanning  the  range  [0  •  l],A^(p) 
is  increasing  for  correlations  less  than  p*  and  decreasing  for 
correlations  greater  than  p*  when  k  >  1.  Ra*  the  lower  Frechet 

distribution  ^Ctp)  »  exp  -  1  du  for  p  ^  (l-Y) 

f  p  u  -  1 

Q  otherwise. 

Table  3  shows  the  value  of  Aj^(p)  x  i00%  for  p  =  ,7,.  .5,  .3  for  the 
two  Frechet  distributions.  For  k  >  1,  the  naximum  value  under 
the  Gijnbel  C  model  is  also  given.  As  in  the  parametric  estimation 
problem  the  largest  errors  are  incurred  vhen  k  <  1.  In  all  cases 
the  effect  of  a  deparUire  from  independence  is  the  largest  v4\en 
p  is  small  (i.e.  for  large  t).  The  effect  decreases  as  k 
increases  reflecting  the  fact  that  vhen  k  Xj^  >  >  X2  the  majority 
of  the  system  failures  are  due  to  tiie  failure  of  the  first 
component. 


5.  Conclusions 

The  results  presented  in  tiiis  paper  diow  that  for  all  five 
bivazoate  exponential  models  one  may  be  apprecialily  misled  by 
falsely  assuning  independence  of  ccn^nent  lifetimes  in  a  series 
system.  The  amount  of  error  incurred  in  modeling  system  reliability 
not  only  depends  upon  the  correlation  be^ieen  component  lifetimes 
but  also  on  the  level  of  system  reliability.  The  error  in 
flodeling  mean  system  life  similarly  depends  upon  the  correlation 
and  the  length  of  maan  system  life.  Both  quantities  depend  on  the 
relative  ma^rtitudes  of  the  parameters. 

For  the  dual  problem  of  estimating  oooipQnent  reliability 
based  on  data  ftxa  a  series  system,  it  appears  that  departures 
from  independence  are  of  greater  oonsequence.  Both  parametric 
and  nonparametrLc  estuoators  of  relevant  ccnponent  parameters 
are  inconsistent.  Bias  increases  dranatically  as  the  correlation 
gets  further  from  zero.  Hoviever,  the  five  models  do  not  exhibit 
appreciable  differences  in  bias  and  mean  sqimured  error  as 
oorrelation  (Ganges.  This  suggests  that  these  models  may  belong 
to  a  large  class  of  bivariate  exponential  distributions  which 
possesses  the  properties  exhibited  here. 
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Appendix 
Gumbel  A 

Gumbel  8 

Gumbel  C 

Downton 


2  -  Foimilas  for  p  s  POCj^  <  X^) 

-  KX^  <  X,)  =  1/2  .  j,/i-  «p(  ^ 

v4iere  «  C  * )  is  the  survival  function  of  a  standard  nomal 
random  variable.  (A2.1) 


-  P(X^  <  X^)  =  Xj^  *  Up 

7xpxp7xpn?xpr2Xj^+x2) 

(A2.2) 

-  P(Xj^  <  X^)  s  X® 

Cxjn®)  (A2.3) 


P(X^ 


X^) 


=  2Xj^X^.(l-p) 

-  **0  X3^X2)(Xj^-X2  ♦  v(xj^n2)^  -  '♦pXj^X2) 

(A2.4) 


Oakes  -  .°(X,  <  X.)  found  numerically. 
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^pendix  3  - 


Gumbel  A  -  Hj^(x)  s  expC-  Xj^x  - 


ajmbel  B  -  Hj^(x)  =  exp  -  [  Cl  +  •♦oCl-expC-A  t))(l-2expC-A,t))]  (A3. 2) 

J  - = - = - dt 

Q  Cl  +  4o(l-exp(-X  t))(l-exp(-Aj^t))  ] 


Gumbel  C  -  Hj^(x)  =  expl-  X™ 


(A3. 3) 


(X™vX^) 


Dawnton  -  Found  numerically  due  to  no  close  farm  solution  for 
F(t). 


Oakes  -  Hj^(x)  =  exp  -  |  Xj^exp(Xj^(e-l)t) 


{exp(Xj^(6-l)t)  *■  expfX^Ce-Dt)  -1} 
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Asymptotic  Modeling  Error  in  Estimating  U'j^  for  K  *  1 
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Asymptotic  Error  in  the  Nonparametric  E^stimator*  of  the 
First  Component  Survival  Function  at 
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FIGURE  6C 

Ar;vmt>tutif:  Error  in  the  Nonparametric  estimator  of  the 
First  Cnm|jonent  Survival  F'inrrii>n  at 

tp  =  .25,  K  =  .67. 
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ABSTRACT 

A  random  environmental  effects  model  is  proposed  for  conq)eting  risks  experiments.  The 
model  assumes  a  random  stress,  Z,  which  changes  the  scale  parameter  of  each  of  the  assumed 
Weibull  times  to  occurrence  of  the  risks.  Some  general  properties  of  the  model  are  discussed,  and 
specific  properties  for  a  Uniform  or  Ganuna  stress  naodel  are  presented.  Estimation  of  parameters 
under  the  Gamma  stress  model  is  considered,  and  a  new  estimator  based  on  the  scaled  total  time  on 
test  transform  is  presented. 

INTRODUCnON 

The  problem  of  competing  risks  arises  naturally  in  a  number  of  engineering  or  biological 
e;q>eriments.  In  such  experiments,  for  some  items  put  on  test,  the  prinuuy  event  of  interests  (such 
as  death,  component  failure,  etc.)  is  not  observable  due  to  the  occurrence  of  some  competing  risk 
of  removal  from  the  study  (such  as  censoring,  failure  from  a  different  component,  etc.). 

Competing  risks  arise  in  an  engineering  context  in  analyzing  data  from 

(a)  series  systems, 

(b)  field  tests  of  equipment  with  a  fixed  test  time  and  a  random  or  staggered  entry  into 
the  study,  or 

(c)  systmns  witii  multiple  failure  modes. 

Conq)eting  risks  arise  in  biological  q)plications  in  analyzing  data  frxnn 
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(a)  cliakal  trials  widi  a  fixed  trial  duration  and  staggered  entry 

(b)  clinical  trials  with  some  patients  withdrawing  fiom  the  trial  prior  to  response 

(c)  studies  of  the  time  to  death  fix>m  a  variety  of  causes 

A  common  assumption  made  in  analyzing  competing  risks  experiments  is  that  the  potential 
(unobservable)  times  to  occurrence  of  the  con^ting  risks  are  independent  This  assunq)tion  is  not 
testable  due  to  the  identifiability  problerrt  That  is,  for  any  dependent  conqjeting  risks  model,  there 
exists  an  independent  competing  risks  model  whkh  yields  the  same  observables.  (See  Basu  and 
Klein  (1982)  for  details.)  However,  Moeschberger  and  Klein  (1984)  show  that  an  investigator  can 
be  appreciably  misled  in  modeling  competing  risks  by  erroneously  assuming  indepencence. 

In  this  paper  we  present  a  model  for  dependence  between  the  various  risks  by  assuming  that 
dependence  is  due  to  some  common  environmental  factor  which  ^ects  the  potential  times  to 
occurrences  of  each  risk.  In  section  2  we  present  die  model  and  shufy  its  pix^)erties  for  bivariate 
series  and  parallel  systems.  In  section  3,  we  consider  estimation  of  the  noodel  parameters  for 
competing  risks  systems. 

2.  THE  MOTEL 

For  simplicity  we  shall  consider  tire  problmn  of  bivariate  systems  and  discuss  our  model  in 
terms  of  engineering  applications.  We  assume  that  under  ideal,  controlled  conditions,  as  one  may 
encounter  in  the  laboratoiy  in  the  testing  or  design  stage  of  development,  the  time  to  failure  of  the 
two  comptments,  to  be  linked  in  a  system,  are  Xq  and  Yq.  We  suppose  that  under  these 

conditions,  Xq,Yq  have  survival  functions  Fq,  Gq  on  [0,  <>•).  We  assume  that  both  Xq  andYQ 

follow  a  WeibuU  form  with  parameters  (T]  j,  Xj)  and  (  T]2>  ^)»  respectively.  That  is,  Fq(x)  =  exp(- 

X^x^l).  The  WeibuU  distribution,  which  may  have  increasing  ( r|  >1),  decreasing  (t^  <  1)  or 

constant  failure  rate  (t] »  1)  has  been  shown  experimentaUy  to  provide  a  reasonable  fit  to  many 
different  types  of  survival  data.  (See  Bain  (1978)).  We  now  link  the  two  components  into  a 
system  in  such  a  way  that  under  ideal  lab  conditions  the  two  components  are  independent 

Now  suppose  that  the  above  system  (Xq,  Yq)  is  put  into  operation  under  usage  conditions. 

We  suiq;>ose  that  under  such  conditions  the  effect  of  the  environment  is  to  degrade  <x  improve  each 
conqxment  by  the  same  random  amount  That  is,  the  effect  of  the  enviixmment  is  to  select  a 
random  factor,  Z,  from  some  distribution,  H,  which  changes  the  maginal  survival  functions  of  the 
two  components  to  Fq^  and  Gq^.  A  value  of  Z  less  than  one  means  that  component  reliabilities  are 
simultaneously  inqnoved,  while  a  value  of  Z  greater  than  one  impUes  a  joint  degradation.  The 
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resulting  joint  reliability  of  the  two  conqxnients'  lifetimes,  (X,Y)  in  the  operating  environment  is 

F(x,y)  =  E[exp(-Z(X|X^l+  )].  (2.1) 

This  model  has  been  proposed  by  Lindley  and  Singpurwalla  (1984)  in  the  reliability  context 
when  Fq,  Gq  are  exponential  and  H(  )  follows  a  gamma  distribution.  This  basic  dependence 

structure  was  also  proposed  by  Qayton  (1978)  to  model  associations  in  bivariate  survival  data,  and 
later  by  Oakes  (1982)  to  model  bivariate  survival  data.  Hutchinson  (1982)  proposed  a  similar 

model  when  H(  )  has  a  gamma  distribution  and  FQ(t) »  GQ(t) »  exp(-tn). 

The  model  described  above  for  a  general  distribution  of  the  environmental  stress  has  a 
particular  dependence  structure  which  we  summarize  in  the  following  lemmas. 

Lemma  1.  Let  (X,Y)  follow  the  model  (2.1)  wl^re  Z  is  a  positive  random  variable  with  finite 

r  s 

( —  +  — ^)“  inverse  moment  Then 
^2 


i 

s 


m 


i 


.  ^ 

'I  . 


S  -f: 


I’l} 


m 


-r/Tij  ■8/n2  "(r/Tii  4.s/n2) 

E(xr  Y*)  *  Xi  r(l  +  r/Hi)  r(l  +  s/7i2)  E(Z  ) 


(2.2) 


The  proof  follows  by  noting  that  given  Z  »  z,  (XY)  are  independent  Weibulls  with  parameters 

,  -r/tli  -Ttr\i 

(7)  1,  Xj  z)  and  (7l2.  X2  z),  respectively  and E(X*'|Z*z)  =  Xj  z  r(l+  rA)  j)  with  a  similar 

expression  for  Y^.  When  the  appropriate  moments  exist  we  have 

(A)  E(X)  =  E(Xo)E(Z  S, 

-lAll  -lAli 

(B)  V(X)  =  E(Xo2)Var(Z  )  +  E(Z  )2var(Xo), 

-lAli  -l/tl7 

(C)  Cov  (X,Y)  =  E(Xq)  E(Yq)  Cov(Z  ,Z  )  which  is  greater  than  0. 


If  11 1  =  il2  =  ^  correlation  between  (X,Y)  is 


r(l+  l/n)^  Var(Z'I^) 


var(z-i^)  r(i+2Ai)+  (r(i+2Ai)  -  r(i+i/n)^)  E(r  i^)2 
In  this  case  the  correlation  is  boundedabove  by  m+lA])^  /  r(l+2Ar|).  Figtue  1  shows  the 
maximal  correlation  as  a  functcion  of  11  for  t)  €  (0, 10).  Note  that  this  maximal  cmrelation  is  an 


increasing  function  of  T].  One  can  also  show  that  F(x,  y)  is  positive  quadrant  dependent  for  any 


Exact  expressions  for  competing  risks  quantities  of  interest  can  be  computed  when  a 
particular  model  is  assumed  for  the  distribution  of  Z.  We  shall  consider  the  gamma  and  uniform 

models.  Consider  first  the  gamma  model  with  h^(z)  =  b®  exp(-bz)/r(a),  z  >  0.  For  this 
model,  the  joint  survival  function  is 


b^ 

F(x,y)  = -  (2.3) 

[b+  X^x  +  r 

which  is  a  bivariate  Burr  Distribution  (see  Takahasi  (1965)),  the  marginal  distributions  are 
univariate  Burr  distributions  with 


E(X)  =  (X^/b)  ^^^(1+1/T\i)r(a-  1/Tii)  r(a),  if  a  >  1^, 

-2Aii  r(i+2Aii)r(a-2A\i)  r(i+i/ni)r(a-i/Tii) 

Var(X)  =  (Xi/b)  { - 1 - L  -[ - i - L]2},ifa>2Aii 

r(a)  r(a) 

with  similar  expressions  for  E(Y),  Var(Y).  The  covariance  of  (X,Y)  is 


-i/Tii  -i/Ti2  na-i/ni-i/n?)  na-i/n?)  na-i/ii^) 

Cov(X,Y)  =  (Xj/b)  (X2/b)  r(l+l^l)r(l+l/Ti2){ - - - } 

r(a)  r(a) 


for  a  >  l/q  2  +  1/Ti2>  For  the  gamma  model,  the  reliability  function  for  a  bivariate  series  system  is 
given  by 

^1  ^2  , 

Rs(t)  =  (l+(Xi/b)t  +(X2/b)t  (2.4) 

and  for  a  parallel  system  by 

Rp(t)  =  (l+(  Xi/b)t^  +  (1+  (X2/b)t^^)'*  -  (l+(  Xi/b)t^^+  (X2/b)t^^)-a  (2.5) 

Hgures  2A-E  arc  plots  of  the  series  system  reliability  for  Xj»  1,  X2  -  2  and  sevo^  combinations 


of  T]  j,  112.  Each  figure  shows  the  reliablity  for  a  =  1/2, 1, 2, 4,  and  the  independent  Weibull 

model.  In  all  cases,  b  =  1.  For  these  figures  we  note  that  for  fixed  ^  'Hi*  'n2>  ^  series 

system  reliability  is  a  decreasing  function  of  the  shape  parameter  a.  Figures  3A-E  are  plots  of  the 
parallel  system  reliability  (2.5)  for  the  above  parameters.  Again,  the  reliability  is  a  decreasing 
function  of  a.  Also  in  both  the  series  and  parallel  system  reliability,  the  sh^  of  the  reliability 
function  is  quite  different  from  that  encountered  under  independence. 

The  gamma  model  is  a  reasonable  model  for  the  environmental  stress  due  to  its  flexibility  and 
the  tractability  of  the  model  in  obtaining  close  form  solutions  for  the  relevant  quantities  and  in 
estimating  parameters.  However,  in  some  cases,  such  as  when  the  operating  environment  is 
always  more  severe  than  the  laboratory  environment,  the  support  of  H  may  be  restricted  to  some 
fixed  interval.  A  possible  model  for  such  an  environmental  stress  is  the  uniform  distribution  over 
[a,b].  For  this  model,  the  joint  survival  function  is 


Til  ^2  ^1  ^2 

[exp  (-b(A.i  X  +  A,2y  ))  -  cxp(-a(  X^x  +  A^y  ))] 
F(x,y)  =  - 

(b-a)(Xix  +  ^2  y  ) 

-1/Tll  (ril-l)/Tll  (ili-l)/tii 

E(X)=Xi  r(l+l/ili)Tii(b  -a  )/{(ili-l)(b-a)] 

=  lln(b/a)/[Xi(b-a)] 


(2.6) 


if  Tl|^  ^  1 
ifTli=l, 


-2/ Til  (Tli-2)/tli  (t1i-2)/i1i 

Var(X)  =  TliXi  {r(l+2/Tli)ili(b  -a  ^ 


(b-a) 


-  (Tll-l/Tli)  (Tli-l/Tli) 

-r(l+l/ni)^Tii(b  -a  )2  if  111  1,2 

(Tli-l)2(b-a) 


2/(Xi2ab)  -  lln(b/  a)2/[(b-a)  Xil^ 


ifTii=l 


^  l!n(b/  a)  n 

^1  ^  TTT""™ '  T75?r 
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FIGURE  2  fl 

SERIES  SYSTEM  RELIABILITY  UNDER  GAMMA  (fl,l)  MODEL 
FOR  THE  ENVIRONMENTAL  STRESS. 

>>1=1.0,  X2=2.0,  ^»i  =  1.0.  112=1.0. 
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FIGURE  2  B 

SERIES  SYSTEM  RELIABILITY  UNDER  GAMMA  (A, 1)  MODEL 
FOR  THE  ENVIRONMENTAL  STRESS. 

>^1  =  1.0,  X2=2.0,  ’iisl.O.  ^2=2.0. 
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FIGURE  2  C 

SERIES  SYSTEM  RELIABILITY  UNDER  GflMMfl(fl.l)  MODEL 
FOR  THE  ENVIRONMENTAL  STRESS. 

^1=1.0*  X2=2.0»  *>1  =  2.0.  ■’>2=2.0. 
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FIGURE  2  D 

SERIES  SYSTEM  RELIABILITY  UNDER  GAMMA  (A, 1)  MODEL 
FOR  THE  ENVIRONMENTAL  STRESS. 

>‘1  =  1.0,  >‘2=2.0,  ’11  =  1/2,  ’12=1/2. 
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FIGURE  3  B 

PARALLEL  SYSTEM  RELIABILITY  UNDER  GAMMA  (A, 1)  MODEL 
FOR  THE  ENVIRQNMENTftL  STRESS. 

Xj  =  1.0,  >2=2.0,  ’fisl.O.  ’»2=2.Q. 


FIGURE  3  D 

PARALLEL  SYSTEM  RELIABILITY  UNDER  GAMMA  (A. 1)  MODEL 
FOR  THE  ENVIRONMENTAL  STRESS. 

Xi=1.0,  X2=2.0,  ^1=1/2.  ^2=1/2. 
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For  this  model,  the  reliability  function  for  a  series  system  is 

111  112  „  'Hi  ,  ^2 

R5(t) « [exp(-b(  Xjt  +X2*  )-cxp[-a(Xjt  +^2!  ) 


issm 


<•* 


.  ,  ^2 
(b-a)(  Xjt  +  ^2 1  ) 


and  for  a  parallel  system  is 

Til  ■ni  ^2  ^2 

Rn(t)  =  [exp(-b(Xi  t  )  -exp(>a  Xi  t  )  +  [exp(-b  Xm  )  -exp(-a  Xm  ) 

^  - - - - —  - = - 1 -  -Rs(t) 

(b-a)  Xjt  (b-a) 

Figures  4A-E  show  the  reliability  for  a  series  system  and  figures  5A-E  for  a  parallel  system 

under  the  uniform  model  for  various  combinations  of  X^,  X2,  “n  i>  ^  2  Notice  that  when  A  = 

.25,  B  3 .75,  which  corresponds  to  an  operating  environment  which  is  less  severe  than  the  test 
environment,  the  system  reliability  is  greater  than  that  expected  under  independence,  while  when 
(a,b)  3  (1.25, 1.75)  or  (1.,  2),  which  corresponds  to  an  environment  more  severe  than  the  test 
environment,  the  system  reliability  is  smaller.  Also  when  the  (a,b)  contains  1,  which  corresponds 
to  an  environment  v^h  incurs  the  possibility  of  no  differential  effect  from  that  found  in  the 
laboratOTy,  there  is  little  difference  in  the  dependrat  and  independent  system  reliability. 

3.  Estimation  of  Parameters  Under  Gamma  Model 

Consider  the  model  (2.3)  with  tij  3  t|2  =  T|.  For  this  model,  the  reliablity  for  a  series  system 
is 

(Xi+Xo) 

Rjj(t)  3(1+ - t^)-a.  (3.1) 

b 

Notice  that  this  model  depends  only  on  two  parameters  0  =  (Xj+  X^lh  and  a  so  that  if  we  had  data 
only  from  systems  (xi  test  in  the  operating  environment,  the  only  identifiable  parameters  are  a,  0, 


FIGURE  4  fl 

SERIES  SYSTEM  RELIABILITY  UNDER  UNIF  (fl, 
FOR  THE  ENVIRONMENTAL  STRESS. 

Xj=1.0,  X2=2.0,  ’»]  =  !. 0,  n2=l.Q. 
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FIGURE  4  B 


SERIES  SYSTEM  RELIABILITY  UNDER  UNIF(A,B)  MODEL 
FOR  THE  ENVIRONMENTAL  STRESS. 

>‘1  =  1.0,  ^2=2.0,  »n  =  l.Q.  ^2=2.0. 

o 

o 


TIME. 
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FIGURE  5  B 

PARALLEL  SYSTEM  RELIABILITY  UNDER  UNIF(R,B)  MODEL 
FOR  THE  ENVIRONMENTAL  STRESS. 

Xi  =  l.O,  X2=2.Q,  ’iisl.O.  n»=2.Q. 


FIGURE  5  D 

PARALLEL  SYSTEM  RELIABILITY  UNDER  UNIF(R,B)  MODEL 
FOR  THE  ENVIRONMENTAL  STRESS. 

Xi=1.0,  X2=2.0,  Vt=l/2,  ^2=1/2. 


FIGURE  5  E 

PARALLEL  SYSTEM  RELIABILITY  UNDER  UNIFIA.B)  MODEL 
FOR  THE  ENVIRONMENTAL  STRESS. 

Xi=1.0,  X2=2.0.  Vi=l/2.  ^2=2.0. 


T],  not  Xj,  X2,  T],  a,  b.  However,  in  many  instances  we  have  extensive  data  on  the  performance  of 


the  conqxments  in  the  lab  under  ideal  operating  conditions  so  that  one  may  ctxisider  X^  X2,  Ti  to  be 

known  based  on  estimates  from  this  data.  We  shall  focus  on  die  problem  of  estimating  6  and  a, 
based  on  data  on  die  system  failure  times  collected  in  the  operating  environment  Let  t^, t^  be 

the  failure  times  for  n  such  systems  put  on  test  and,  let  w^  =  i  =  1, ...,  n. 

Prior  to  attempting  to  estimate  (a,  6),  we  would  like  to  check  if  the  model  (3.1)  is  feasible. 

A  graphical  check  of  this  model  can  be  done  through  die  scaled  total  time  on  test  (STIOT)  plot  of 
Barlow  and  Campo  (1975).  The  STTOT  fw  W  is 

,F-l(t) 

J  Rs(t)dt 

Gw(t) - , - *  l-(l-t)(^^-lV«  for  a  >  1.  (3.2) 

f1(1) 

/  Rs(t)dt 

0 

Note  that  (3.2)  depends  only  on  a.  Figure  6  shows  the  form  of  the  STTOT  for  several  values  of  a. 

Notice  that  for  all  a,  the  STTOT  is  below  the  45”  line  (which  corresponds  to  exponential  system 
life)  since  the  hazard  rate  of  the  series  system  is  decreasing.  Let 
i 

Tn(W(i))  -  2  W(j)  +  (n.i)W(i),  (3.3) 

j=l 

where  Wq^  ^  ^(2)  ^  •••  ^  ^(n)  ordered  systems  failure  times  be  the  total  time  on  test  at 

The  empirical  STTOT  plot  then  plots  (i/n,  T„(W^))/T(W(„)))  which  can  be  compared  to 
figure  6  for  a  graphical  check  of  the  model  Also,  crude  estimates  of  a  can  be  obtained  by 
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(X>inparing  the  empirical  and  theoretical  STTOT  plots.  When  there  is  no  random  environmental 
effect  and  the  conq)onents  are  independent,  then  the  empirical  STTOT  plot  should  look  like  the  45° 


line.  Also  as  a  tends  to  infinity  this  plot  approaches  the  45**  line. 


We  now  consider  several  estimates  of  a  and  9.  The  log  likelihood  for  the  model  (3.1),  based 
on  a  sample  of  size  n,  is 

n 

L(a,6)  =  nltna  +  niln9  -  (a+l)Z  in  (1+0  W:)  (3.4) 

i=l 


so  that 


3/9a  L(a,6)  =»  n/a  -  L  in  (1  +  0  Wj) 
i=l 


and  9/30  L(a,O)  =  n/0-(a+l)  Z  wj/(l+0w:) 

i=l 


For  (3.5)  we  note  that  the  maximum  likelihood  estimator  of  a 


%nle“ 


I  in  (1  +  0  Wj) 
i=l 


and  the  maximum  likelihood  estimator  of  9  is  the  solution  to 


n  n  Wj 

_  -( - +1)(2 - )-0. 

0  Zin(l+0Wi)  l  +  0Wi 


One  can  show  that  0  is  positive  if  n  Z  w;^  >  2(  Z  w:)^. 

i-1  i»l 


In  such  case  9jQ]g  is  obtained  by  solving 


(3.8)  numencally. 


A  second  estimator  of  (a,  0)  is  the  mediod  of  moments  (mme).  Since  E(W)  =  [  6(a-l)]' 


and  E(w2)  *  2[  e2(a-l)(a-2)]-l  where  a  >  2,  we  have 


Swi2 

S  -  2(2 


^  -= — ^ ®nime 


1  Wi2.2(2;wi)2 

2  wjd  Wi)2 


(3.11) 


provided  that  (3.9)  holds.  If  (3.9)  does  not  hold,  tten  this  estimator  does  not  exist 

A  diird  estimator  was  suggested  by  Berger  (1983)  in  a  different  context  Hie  suggested 


estimating  6  a  modified  methods  of  moments  estimator  6|^  »  (aw)-l.  (3.12) 


where  w  =  2  wj/n. 


which  is  used  as  the  true  value  of  9  in  the  likelihood  (3.4)  so  tibat  the  estimator  of  a  is  the  solution 


A  final  estimator  is  based  on  the  STTOT  plot  Let  »  fln(l-i/n)  and  : 

lln(l-Tn(W(jp/Tn(W(jjp),  i  =  1, ...,  n-1.  If  (3.2)  holds,  then  we  should  have 
Dj » (l-l/a)Ci,  i  =  1, ...,  n-1. 


so  the  value  of  a  which  minimizes 


n-1  - 

2  (D{  -  (1-1/a)  C{r  is  a  reasonable  estimator  of  a 
i»l 


The  resulting  estimator  is  a^g  ^ 


2q2-2qDi 


which  is  in  the  parameter  space  if  Z  >  Z  Ci  Dj.  A  better  estimator  should  be  obtained  by 
weighting  the  D^'s  differently  since  for  i  <  j,  Var  (D^)  <  Var  (Dj).  The  variance  of  Dj  depends  on 
the  unknown  parameter  a  so  we  weight  by  the  variance  of  computed  under  an  assumed 
exponential  distribution.  The  variance  of  in  that  case  is 
i  1 


V|  —  Z  ■■■  "  A ' » 1  ^  1»  •••»  n*! 
j=l  (n^ 


(3.16) 


so  that  die  weighted  least  squares  estimator  of  a  is 


^vls  = 


Iq2/Vi 


(2—  -2 


qq 


.ifZq2/Vi>ZCiDiA^i. 


(3.17) 


(Dnce  we  have  obtained  a  by  either  of  the  two  least  squares  estimators,  we  substitute  this  value  into 
(3.6)  and  solve  this  equation  numerically  for  6^  or 

The  condition  Z  >  Z  q  q/Vf  includes  a  few  more  possible  sanqiles  than  the 
condition  (3.9)  for  the  other  three  estimators.  However,  those  samples  which  satisfy  Z  q^/V^  > 
Z  q  /V^  for  which  (3.9)  fails  to  obtain  yield  very  large  estimates  of  9.  Since  a  reasonable 
model  for  T  when  6  and  a  are  not  estimable  is  the  independent  Weibull  series  system  which  has 
system  reliability  very  close  to  (3.1)  when  a  is  very  large,  this  is  not  a  problem.  Figures  7a  and  7b 
are  scaled  total  time  on  test  plots  from  two  simuilated  samples  of  size  30  from  (3.1)  with  a  -  3, 9  = 
1.  Looking  at  figure  7a,  we  see  that  the  estimated  scaled  total  on  test  doesn't  look  too  different 
fiom  the  45**  so  that  an  exponential  model  might  not  be  unreasonable.  For  this  data  set  only  the 
weighted  least  squares  estimator  exists  and  it  yields  a^;^x.5  »  45.33  and  9  » .0567.  For  the  data  in 
figure  7b  all  estimates  exist,  and  we  have 
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‘  C 


®mle  “  ®inle  “ 

^mme  *  a„,„e  .  4.86 

®ber  “  ®ber  “  7.02 

01s  =  .739  ais  -  3.58 

®wls  “  *^70  ^Is  “ 

To  study  the  properties  of  these  estimators,  a  small  scale  Monte  Colo  study  was  performed. 
Random  samples  of  size  n  =  15, 30, 50, 75,  or  100  were  generated  with  Xi  +  X2  =  3,  b  =  3,  so  0 
=  1  and  a  a  2, 3, 5.  1000  samples  were  generated  for  each  combination  of  n  and  a.  The  bias, 
standard  deviation  of  the  estimates  and  n,  the  number  of  samples  where  the  estimator  exists  is 
reported  in  table  1  for  a,  table  2  for  0,  and  in  table  3  for  an  estimator  of  the  system  reliability 
obtained  from  (3. 1)  at  =  9.085.  The  true  system  reliability  at  tg  is  .8255  when  a  =  2,  .75  when  a 
»  3,  and  .619  when  a  =  5.  Also  reported  in  each  table  is  the  bias  and  standard  deviation  of  the  least 
square  and  weighted  least  square  estimators  when  tiiey  are  restricted  to  those  samples  where  the 


^mme  * 


other  estimators  exist 


From  these  tables  we  note  that  Berger's  modified  estimatcr  performs  very  poorly.  Also  the 
weighted  least  squares  estimator  allows  for  estimation  of  parameters  in  many  mcne  samples  when  n 
is  small.  In  general  the  maximum  likelihood  estimator  outyerforms  the  other  estimators,  however, 
when  the  weighted  least  squares  estimator  is  restricted  to  tiiose  sanq>les  where  the  maximum 
likelihood  estimator  exists,  this  estimator  performs  much  better  when  n  is  stnalL  The  somewhat 
better  performance  of  the  MLE  in  terms  of  Itias  is  deceptive  since  some  of  the  estimates  of  a  are  less 
than  one,  which  implies  tiiat  the  mean  system  reliability  is  infinite.  Also  the  weighted  least  squares 
estimator  of  system  reliability  seems  to  ouq)erf(xm  the  (AhCT  estimators  of  the  system  reliability  in 
spite  of  its  relatively  poor  performance  as  an  estimator  of  0.  Our  recommendation  is  to  use  the 
weighted  least  squares  estimator  since  it  more  often  provides  estimators  of  tiie  relevant  parameters 
and  is  somewhat  easier  to  compute. 
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0. 

511 

978 

-.345 

0. 

297 

978 

-.716 

0. 

104 

978 

055 

0. 

,465 

952 

- 

.161 

0. 

509 

522 

1. 

352 

3. 

109 

665 

0. 

715 

1. 

,609 

558 

0 

.578 

1. 

536 

522 

-.238 

0. 

601 

522 

-.546 

0. 

348 

516 

1. 

084 

1 . 

,599 

458 

0 

.827 

1. 

561 

674 

0. 

558 

1. 

609 

752 

0. 

366 

1. 

.118 

669 

0 

.  180 

1. 

026 

674 

-.199 

0. 

576 

674 

-.584 

0. 

282 

660 

0. 

523 

1 . 

,104 

601 

0 

.286 

1. 

029 

801 

0. 

256 

1 . 

559 

850 

0. 

184 

0, 

,869 

756 

0 

.  105 

0 . 

863 

801  ■ 

-.  1'’3' 

0. 

541_ 

801 

-.615 

I>. 

232 

787 

0. 

267 

0. 

850 

722 

0 

.150 

0. 

358 

GO 

0. 

129 

0. 

917 

915 

0. 

112 

0, 

,728 

327 

0 

.014 

0. 

747 

893 

-.189 

0. 

535 

893 

-.628 

0. 

21'.' 

878 

0. 

153 

0 . 

.715 

614 

0 

.020 

0. 

744 

992 

0. 

033 

0. 

683 

913 

0. 

020 

0, 

.623 

8.:.3 

- 

.  064 

0. 

666 

892 

-.206 

0. 

494 

892 

-.644 

185 

879 

0 . 

055 

0. 

.615 

821 

- 

.050 

0. 

663 

TABLE  3 

BIAS  AND  STANDARD  DEVlATiaN(SO»  QF  ESTIMATORS  OF  SYSTEM  RELIABILITY  AT 
MAXIMUM  WEIGHTED  METHOD  OF 

LIKELIHOOD  LEAST  SQUARES  LEAST  SQUARES  MOMENTS 


T».9085 

BERGER’S  METHOD 


A 

N 

M 

BIAS  SO 

M 

BIAS 

SD 

M 

BIAS  SD 

M 

BIAS  SD 

M 

BIAS 

SD 

2 

■■'is 

'769 

-.012  .0647 

'832 

-.004 

.0386 

‘762 

0.002  .0388 

'770 

0.037  .0503 

"770 

0.064 

.0463 

766 

-.006 

.0577 

715 

0.002  .0386 

2 

30 

916 

-.005  .0473 

953 

0.002 

.0424 

877 

0.010  .0434 

916 

0.037  .0357 

916 

0.069 

*  0  ^  5 

912 

0.001 

.0426 

837 

0.009  .0434 

O 

50 

979 

-.001  .0372 

989 

0.003 

.0349 

936 

O.OlO  .0359 

979 

0.035  .0300 

979 

0.071 

.0233 

976 

0.003 

.0348 

952 

0.010  .0359 

2 

75 

996 

0.000  .0290 

998 

0.004 

.0274 

974 

0.010  .0292 

996 

0.034  .0244 

996 

0.0/2 

.0238 

996 

0.003 

.0274 

972 

0.010  .0293 

2 

100 

999 

O.OOl  .0243 

1000 

0.004 

.0234 

989 

0.010  .0248 

999 

0.034  .0223 

999 

0.0  75 

.0216 

999 

0.004 

.0233 

989 

O.OlO  .0248 

3 

15 

642 

-.018  .0813 

753 

-.010 

.0767 

653 

-.a06  .0748 

643 

0.031  .0661 

643 

0.06  7 

.0616 

636 

-.015 

.0764 

373 

-.008  .0751 

3 

30 

809 

-.007  .0377 

870 

-.003 

.0532 

769 

0.002  .0531 

810 

0.024  .0490 

809 

0.062 

.0472 

804 

-.006 

.0544 

731 

0.001  .0550 

3 

30 

916 

-.003  .0429 

933 

-.001 

.0412 

864 

0.005  .0431 

916 

0.022  .0366 

916 

0.066 

.0340 

912 

-.002 

.0411 

831 

0.003  .0432 

3 

73 

963 

-.001  .0372 

977 

0.000 

.0336 

925 

0.007  .0374 

963 

0.021  .0327 

963 

0.067 

.0313 

958 

0.000 

.0357 

923 

0.007  .0375 

3 

too 

9788 

O.OOl  .0309 

989 

0.002 

.0299 

956 

0.008  .0311 

978 

0.019  .0267 

978 

0.067 

.0261 

978 

0.002 

.0299 

952 

0.008.  .0312 

5 

15 

52.0 

-.0’2.9  .1011 

665 

-.024 

.0967 

558 

-.020  .0977 

522 

0.012  .0926 

522 

0.054 

.0921 

516 

- .  030 

.0968 

438 

-.022  .0968 

5- 

30 

674 

-.022  .0691 

732 

-.020 

.0674 

669 

-.013  .0651 

674 

O.OOl  .0613 

674 

0.045 

.0608 

660 

-.027 

.0665 

601 

-.016  .0652 

3 

50 

801 

-.006  .0545 

830 

-.006 

.0330 

736 

-.002  .0532 

801 

O.OlO  .0494 

801 

0.055 

.0485 

787 

-.  008 

.0528 

722 

-.002  .0535 

3 

75 

893 

-.005  .0442 

915 

-.005 

.0436 

827 

-.002  .0431 

893 

0.007  .0406 

893 

0.051 

.0380 

B78 

-.006 

.  0430 

814 

-.002  .0432 

ir 

%j 

100 

892 

-.002  .0375 

913 

■•.001 

.0392 

835 

0.003  .0380 

892 

0.008  .0353 

892 

C.052 

.  0345 

879 

-.002 

.  0390 

821 

0.003  .0381 
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1.  Introduction 

The  problea  of  coapeting  risks  arises  naturally  in  a  nuaber  of  contexts, 
naaely  Che  aodeling  of  series  systcas  in  reliability,  Che  problea  of  esclaacion 
with  censored  data  and  Che  analysis  of  physical  or  biological  systeas  with 
multiple  failure  modes.  A  cnaann,  uncestable  assuaption  is  usually  aade  that 
Che  potential  failure  Claes  for  each  risk  are  statistically  independent  (see 
Basu  and  Klein  (1982)).  Moeschberger  and  Klein  (1984)  show  chat  an  investiga¬ 
tion  may  be  appreciably  aisled  in  aodeling  series  systeu  reliability  and  in 
estimating  component  paraaeters  by  incorrectly  assuaing  independent  component 
lifetimes.  In  this  paper  we  aodel  dependence  between  coaponencs  through  a 
common  environaental  effect  on  each  component.  Such  a  dependence  structure 
has  also  been  suggested  by  Oakes  (1982).  Lindleyand  Singpurwalla  (198S)  ,  and 
Hutchinson  (1981). 

2.  The  Model 

Consider  a  two-coaponent  system.  Suppose  that  under  ideal,  controlled 
conditions,  as  encountered  in  the  casting  stage,  the  times  to  failure  of  the 
two  coaponencs  ace  X.  and  Y.,  Under  these  conditions,  ***  Independent 

with  marginal  survival  functions,  F.  and  Gq.  Now  suppose  the  two  coaponencs 
are  linked  into  a  syscea  and  exposed  to  Che  environment.  The  effect  of  Che 
environment  is  Co  select  a  randoa  factor,  Z,  from  a  distribution,  H(z)  ,  which 

Z  Z 

changes  the  marginal  survival  functions  of  Che  two  components  to  > 

respectively.  A  value  of  Z  less  Chao  one  aeans  that  component  reliability  is 
improved,  while  a  value  greater  Chan  one  iaplies  a  Joint  degradation.  In  the 
sequel  we  assume  chat  end  Y,,  follow  a  Weibull  distribution  with  parameters 

(a^.A^).  and  (a^^.A^).  That  is.  Fg(x)  -  exp(-Xj^x^).  The  resulting  joint 

reliability  of  the  two  components'  lifeciaes,  (X,Y),  in  the  operating  environ- 

®Y 

menc  is  F(x,y)  ■  E(exp(-Z(AjjX  *  Ay,y  Jf.  F(x,y)  is  positive  quadrant  dependent. 
Also  E(X)  -  E(Xg;E(Z"^^)  ;  V(X)  -  E(Xg2)E(Z'^^®«)  -(E  (Xg)E(2‘^'®*)) ^  ;  and 

Cov(X,Y)  •  E(XQ)E(Yg)  Cov(z"^^‘^,  z"^^®^) .  The  correlation  ia  always  positive 

and  is  bounded  above  by  Fd  l/a)^/r(l  ♦  2/a)  when  «,j  •  «„  "  «.  Explicit, 
chough  lengthy,  formula  for  the  aoaencs  of  (X.Y).  syscea  and  coaponenc 
reliability,  and  syscea  coaponencs  can  be  obtained  when  Z  is  aasuaed  to  be 
either  a  uniform  or  gaasiar  random  variables.  The  gaaaa  case  leads  to 
bivariate  Burr  distributions. 

3.  Estimation  ^ 

The  esc isa cion  of  aodel  parameters  is  carried  out  under  the  assumption 

chat  Z  has  a  gaaaa  distribution  with  density  h(x)as*“*^  exp(-bs).  Since  Che 
parameters  Xy  are  not  identifiable  when  only  data  from  either  a  series  or 


parallel  syseca  is  available,  we  incorporate  sample  information  obtained 
independently  on  each  component  under  test  conditions.  Maximum  likelihood 
and  method  of  momenta  estimators  are  obtained  and  their  properties  are  studied 
by  Monte-Carlo  methods  since  no  closed  form  maximum  likelihood  estimates  are 
available. 
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Resume 


Un  modele  pour  les  systaes,  dependants  dans  1 'analyse  de  la  fiabilite  est 
examine.  Le  modele  suppose  que,  sous  les  conditions  ideales,  les  temps  de 
survie  des  constituents  du  systne  ont  dss  distributions  Ueibull  independantes. 
Sous  des  conditions  d 'operation  un  facteur  exterieur  alcatoire  affecte  cheque 
conatituant  simultandsMnt  an  multipliant  son  taux  de  hasard  par  une  quant  itd 
aleatoire.  Las  proprictes  de  ce  modele  et  I'estlmation  des  parsmetres  du 
modble  sont  considdrds,  k  partir  des  exeaq>les  concrets  du  laboratoire  at  de 
la  pratique. 


